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NOMENCLATURE 


Blowing  parameters 

Stanton  number  for  mass  transfer 

Enthalpy 

Solid  (graphite  or  char  layer)  enthalpy 
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Streamline  metric  coefficient 
Mass  flux 

Thermochemical  graphitic  blowing  rate 

Phenolic  resin  gas  blowing  rate 

Mach  number 

Pressure 

Prandtl  number 

Heat  flux 

Convective  heat  flux 
Re-radiative  heat  flux 
Radial  coordinate 
Reynolds  number 
Nose  radius 

Stream  length  measured  from  the  stagnation  point 
Temperature 
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component  (3IS) 
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Body  normal  velocity  component  (3DMEIT);  circumferential 
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Axial  velocity  component  (3IS) 

Distance  in  surface  normal  direction 
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Y  Polar  coordinate  (2IT) 

Z  Axial  coordinate  (3IS) 

a  Angle  of  attack 

0  Boundary  layer  momentum  thickness;  circumferential  coordinate  (3IS) 

Otj  Body  angle  relative  to  the  centerline 

O'  Circumferential  coordinate  ( 2 I T ) 

p  Viscosity 

p  Density 

4)  Boundary  layer  energy  thickness 
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e  Boundary  layer  edge 

1  Invlscid  flow 

l  Laminar  flow 

r  Recovery  state 

t  Turbulent  flow  or  stagnation  conditions 

te  Stagnation  conditions  at  the  boundary  layer  edge 

tiw  Inviscid  stagnation  conditions  at  the  wall 

tr  Transition  point 

w  Wall 
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SECTION  1 


INTRODUCTION 

For  most  reentry  vehicles  and  missions,  the  aerothermal  environ¬ 
ment  of  primary  concern  for  vehicle  design  is  the  lower  altitude  flight 
regime  (h  <  120  KFT) ,  where  peak  aerodynamic  loads  and  heating  are  ex¬ 
perienced.  Successful  design  and  analysis  of  maneuvering  reentry  vehicles 
require  the  ability  to  accurately  predict  this  aerothermal  environment, 
where  the  Reynolds  number  is  sufficiently  large  that  the  shock  layer  flow 
is  primarily  inviscid,  with  viscous  effects  confined  to  the  thin  boundary 
layer  adjacent  to  the  vehicle  surface  and  to  leeside  regions  of  separated 
flow  at  high  angles  of  attack. 

Sophisticated  flow  field  prediction  techniques  have  been 
developed  that  directly  compute  the  entire  shock  layer,  both  inviscid  and 
viscous  regions,  eliminating  the  need  for  explicit  coupling  between  the 
inviscid  and  viscous  portions  of  the  flow.  These  techniques,  the  para¬ 
bolized  Navier-Stokes  (PNS)  procedures,  are  thus  capable  of  rigorously 
treating  higher  altitudes  (where  the  boundary  layer  is  not  thin),  as  well 
as  leeside  separated  flow  regions.  (A  review  of  the  current  PNS  methods 
was  performed  as  part  of  the  MAT  program  and  is  documented  in  Reference  1.) 

These  detailed  PNS  solutions,  however,  generally  require  a 
large  amount  of  computer  time  and  presently  do  not  constitute  a  design 
technique.  Consequently,  more  efficient  inviscid  flow  field  techniques 
are  often  used  to  predict  the  aerodynamic  forces  and  moments  (e.g.,  normal 
force,  pitching  moment,  and  inviscid  contribution  to  the  axial  force). 
Boundary  layer  calculations  are  used  to  determine  both  the  heat  transfer 


to  the  body  and  the  viscous  shear  forces,  with  the  required  boundary  layer 
edge  conditions  being  obtained  from  the  inviscid  solution.  This  coupled 
inviscid-boundary  layer  solution  procedure  has  been  widely  used  to  obtain 
accurate  flow  field  predictions  for  both  ballistic  and  maneuvering  re¬ 
entry  vehicles. 

For  engineering  calculations,  the  most  appropriate  approach  for 
boundary  layer  predictions  is  the  technique  in  which  integral  boundary 
layer  equations  are  obtained  by  integrating  the  fundamental  boundary 
layer  equations  across  the  layer.  The  solution  of  the  resulting  equations 
determines  the  viscous  parameters  of  primary  interest  -  shear  stress,  heat 
transfer,  boundary  layer  thicknesses,  etc.  -  but  does  not  provide  definition 
of  specific  details  such  as  velocity  and  temperature  profiles  across  the 
viscous  layer.  For  three-dimensional  problems,  the  integral  boundary  layer 
method  is  usually  applied  by  integrating  the  axisymmetric  integral  equa¬ 
tions  along  inviscid  surface  streamlines.  Although  this  approach,  termed 
the  "small  cross-flow  approximation,"  neglects  the  cross-flow  profiles  in 
the  boundary  layer  at  angle  of  attack,  it  has  been  successfully  applied 
to  many  RV  configurations,  providing  reasonably  accurate  boundary  layer 
solutions  with  great  efficiency. 

As  part  of  the  BMO  Maneuvering  Aerothermal  Technology  (MAT)  pro¬ 
gram,  an  existing  state-of-the-art  three-dimensional  inviscid  flow  field 
code  has  been  coupled  to  an  extension  of  an  existing  integral  boundary 
layer  solution  procedure  in  order  to  provide  the  capability  for  efficient 
inviscid/viscous  flow  field  solutions  on  maneuvering  reentry  vehiclts, 
while  taking  maximum  advantage  of  the  inherent  capabilities  of  the  codes 
being  coupled.  The  inviscid  code  used  in  this  coupling  procedure  is  the 


BMO/3 IS  code,  which  is  described  in  References  2-4.  (The  review  of  in- 
viscid  flow  field  codes  for  application  to  MaRV  configurations  conducted 
under  the  MAT  program  is  documented  in  Reference  5.)  The  integral  boundary 
layer  code  used  is  an  extension  of  the  Momentum- Energy  Integral  Technique 
(MEIT),  described  in  References  6-10,  which,  unlike  many  integral  boundary 
layer  techniques,  solves  both  the  momentum  and  energy  integral  equations 
rather  than  invoking  the  Reynolds  analogy  to  relate  heat  transfer  and  skin 
friction. 

The  approach  taken  in  this  coupling  effort  is  described  in 
Section  2  of  this  report,  which  also  details  the  assumptions  made.  The 
inviscid  flow  field  code  used  in  this  effort  is  briefly  described  in 
Section  3,  as  is  the  procedure  developed  in  this  effort  to  compute  the 
inviscid  surface  streamlines.  Details  on  the  3DMEIT  Integral  boundary 
layer  solution  technique  are  provided  in  Section  4. 


SECTION  2 


APPROACH  AND  ASSUMPTIONS 

The  goal  of  this  effort  is  the  development  of  an  efficient, 
robust,  reliable  procedure  for  computing  complete  flow  fields  on  maneuver¬ 
ing  reentry  vehicles  in  the  high  Reynolds  number  flight  regime,  where  the 
thin  boundary  layer  assumption  is  valid.  In  this  regime,  the  Inviscid 
portion  of  the  shock  layer  can  be  accurately  computed  without  regard  to 
the  boundary  layer,  and  the  boundary  layer  solution  can  be  determined 
using  edge  conditions  obtained  from  the  inviscid  solution. 

In  this  task,  the  BM0/3IS  three-o miensi onal  flow  field  code  is 
used  to  define  the  inviscid  shock  layer,  and  an  extension  to  the  MEIT 
integral  boundary  layer  procedure  is  used  to  determine  the  viscous  flow. 
This  extension,  called  3DMEIT,  solves  the  integral  boundary  layer  equations 
along  Inviscid  surface  streamlines  to  efficiently  model  the  three-dimen¬ 
sional  boundary  layer  using  the  "small  cross-flow"  assumption.  The 
coupling  procedure  developed  in  this  effort  provides  definition  of  the 
inviscid  surface  streamlines  and  flow  field  properties  obtained  from  the 
inviscid  solution  in  a  format  compatible  with  the  3DMEIT  boundary  layer 
procedure. 

A  unique  feature  of  this  coupling  procedure  is  the  method  used 
to  define  the  inviscid  surface  streamlines.  In  previous  3-D  Integral 
boundary  layer  solutions,  such  as  those  described  in  References  11  and  12, 
the  Inviscid  streamlines  were  determined  in  an  approximate  manner  using 
the  surface  pressure  distribution  obtained  from  the  inviscid  flow  field 
solution.  In  the  current  effort,  the  inviscid  surface  streamlines  are 
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computed  explicitly  as  part  of  the  invlscid  flow  field  procedure;  this 
technique  is  described  in  Section  3. 

Another  unique  feature  of  the  coupled  3IS/3DMEIT  flow  field 
procedure  is  the  method  used  for  defining  boundary  layer  edge  properties. 
Rather  than  determining  the  edge  entropy  from  a  stream  tube  entrainment 
procedure  (which  works  well  for  axisymmetric  configurations  at  zero  angle 
of  attack,  but  is  difficult  to  apply  in  three  dimensional  problems),  the 
edge  conditions  are  determined  using  a  mass  balance  between  the  invlscid 
and  boundary  layer  flows,  using  the  detailed  information  available  from 
the  finite  difference  inviscid  solution.  This  procedure  is  described  in 
Section  4.2.6. 

A  buffer  code  has  been  developed  to  handle  the  coupling  between 
the  3IS  and  3DMEIT  codes.  This  buffer  code  accepts  as  input  the  3IS 
solution  output  file  and  defines  the  inviscid  surface  streamlines  to  be 
used  by  3DMEIT  in  the  boundary  layer  calculations.  The  format  of  the 
buffer  code  output  is  described  in  Section  3.4. 

Several  requirements  were  established  for  the  3DMEIT  code  in 
this  effort.  These  requirements  were: 

•  the  ability  to  treat  equilibrium  air  thermodynamics, 
in  a  manner  consistent  with  the  3IS  inviscid  code 

d  the  ability  to  automatically  determine  the  wall  tempera¬ 
ture  and  mass  addition  rates  for  ablating  heatshields 
using  steady  state  ablation  models  (or,  as  an  option, 
allow  values  of  Tw  and  m  to  be  input  by  the  user) 

•  the  ability  to  treat  laminar,  transitional,  or  fully 
turbulent  flows. 

These  capabilities  of  the  3DMEIT  code,  and  others  that  have  been  implemented 
are  described  in  Section  4  of  this  report. 


The  coupled  3IS/3DMEIT  procedure  is  applicable  to  any  vehicle 
geometry  which  the  inviscid  (3IS)  code  is  capable  of  treating.  3IS  is 
formulated  in  a  cylindrical  coordinate  system  and  allows  the  definition 
of  generated  afterbody  cross-sections.  Special  geometry  procedures  are 
included  to  allow  simple  specification  of  multi-conics  with  slices, 
cuts  and  flaps.  (Short-comings  do  exist  in  the  current  ability  of  3IS 
to  execute  calculations  over  moderate  to  large  flap  deflection  angles; 
however,  activities  are  now  underway  to  eliminate  this  limitation.) 

The  ability  of  the  3IS/3DMEIT  technique  to  treat  control  surfaces  is 
partially  demonstrated  in  the  sample  calculation  in  Section  4.4.3 
which  involves  a  fixed  yaw  stabilizer  (FYS). 

The  intended  emphasis  of  this  effort  was  to  develop  and 
demonstrate  a  valid  3IS/3DMEIT  coupling  for  afterbody  applications. 
Therefore  it  was  decided  to  Initially  restrict  attention  to  spherical 
nosetip  shapes.  The  axisymmetric  2IT  nosetip  code  is  currently  used 
to  determine  the  inviscid  flow  profiles  associated  with  the  streamlines 
on  a  sphere.  This  limitation  can  be  removed  by  using  the  CM3DT  blunt 
body  code  for  the  treatment  of  non-spherical ,  asymmetric  nosetip  shapes. 
It  is  necessary  to  simply  calculate  the  streamlines  paths  from  a  con¬ 
verged  CM3DT  solution  and  transmit  the  corresponding  inviscid  flow 
field  information  to  the  buffer  code. 

In  addition,  the  coupled  3IS/3DMEIT  flow  field  procedure 
currently  cannot  treat  non-zero  sideslip  angles  because  the  3DMEIT 
geometry  routine  is  restricted  to  the  description  of  bodies  over  a 
half-plane.  The  inviscid  codes  (blunt  body  and  afterbody)  which  are 
used  to  provide  the  inviscid  flow  field  information  necessary  to  drive 
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SECTION  3 


INVISCID  FLOW  FIELD  SOLUTION 

The  1 nvi sc i d  frustum  flow  field  technique  used  in  this  effort  is 
the  BMO/3 IS  code.  Since  spherical  nosetips  are  assumed  in  this  initial 
effort,  the  required  initial  data  for  the  3IS  afterbody  calculations  are 
obtained  from  the  axisymmetric  2IT  nosetip  flow  field  code.  Both  the  2 IT 
and  3IS  codes  are  documented  in  References  2-4  and  are  briefly  described 
in  Sections  3.1  and  3.2  below. 

Sections  3.3  and  3.4  describe  the  buffer  code  developed  to 
handle  the  interface  between  the  1 nvi sc  Id  codes  and  the  3DMEIT  boundary 
layer  procedure  and  the  user  instructions  for  this  buffer  code. 

3.1  BLUNT  BODY  CODE  DESCRIPTION 

Inviscid  flow  field  solutions  on  spherical  nosetips  are  readily 
obtained  with  the  time-dependent  axisymmetric  2IT  code.  This  code  obtains 
the  steady  state  solution  as  the  asymptotic  limit  of  an  unsteady  flow, 
starting  from  an  assumed  initial  flow  field.  The  time-dependent  inviscid 
flow  equations  are  solved  using  the  MacCormack  finite  difference  procedure. 

Formulated  in  a  spherical  coordinate  system,  the  2IT  code  is 
ideally  suited  to  the  calculation  of  flows  on  spherical  nosetips,  since 
the  grid  is  perfectly  aligned  with  the  body  geometry.  Calculations  on 
spheres  with  2IT  are  always  performed  In  wind-fixed  coordinates;  appro¬ 
priate  Initial  data  for  the  3IS  afterbody  code  at  angle  of  attack  are 
determined  by  interpolation  and  rotation  of  the  2IT  solution.  The  2IT 
code  can  treat  either  ideal  gas  or  equilibrium  air  thermodynamics. 
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3.2  AFTERBODY  CODE  DESCRIPTION 


Inviscid  afterbody  flow  field  solutions  are  obtained  for  this 
effort  with  the  3IS  code.  This  code  is  a  forward  marching,  steady, 
finite  difference  flow  field  solution  procedure,  and  requires  only  that 
the  axial  component  of  the  flow  in  the  frustum  shock  layer  be  supersonic 
at  all  points. 

3IS  is  formulated  in  a  cylindrical  coordinate  system  and  allows 
the  definition  of  generalized  afterbody  cross-sections.  Alternatively, 
special  geometry  procedures  have  been  implemented  In  the  code  to  allow 
simple  specification  of  multi-conics  with  slices,  cuts,  and  flaps.  Co¬ 
ordinate  stretching  can  be  used  in  the  3IS  code  to  circumferentially 
concentrate  mesh  points  in  either  the  pitch  or  yaw  planes.  In  addition, 
recent  modifications  to  this  code13  also  permit  radial  clustering  of 
mesh  points  near  the  body  to  Improve  resolution  of  the  flow  gradients 
near  the  wall.  Either  Ideal  gas  or  equilibrium  air  thermodynamics  can 
be  used  In  the  3IS  calculation  procedure. 

3.3  INVISCID  SURFACE  STREAMLINE  CALCULATION 

Inviscid  surface  streamlines  are  defined  for  the  coupled  31 S/ 

3DMEIT  procedure  through  explicit  streamline  calculations  In  the  3IS 

f: 

^  inviscid  afterbody  code.  The  inviscid  streamlines  on  the  spherical  nose- 

t;  . 

;  tip  are  readily  defined  by  noting  that  the  streamlines  remain  in  meridional 

:  planes  in  wind-fixed  coordinates  on  the  sphere. 

It  is  critical  in  this  coupled  Invlscid-boind  ary  layer  approach 
to  ensure  that  the  selected  Inviscid  streamlines  along  which  boundary 
layer  solutions  are  to  be  obtained  provide  complete  coverage  of  the 
••(.  body.  To  accomplish  this  goal,  the  selected  streamlines  Include  those 
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streamlines  that  are  equally  spaced  circumferentially  at  the  initial 
afterbody  plane  (at  the  nosetip  tangency  point)  as  well  as  those  that 
are  equally  spaced  circumferentially  at  the  end  of  the  body.  This  choice 
allows  complete  coverage  of  the  body,  since  at  angle  of  attack  the  equally 
spaced  streamlines  at  the  initial  plane  will  be  swept  around  the  body  and 
will  be  concentrated  near  the  lee  plane  at  the  end  of  the  body,  while 
those  streamlines  equally  spaced  at  the  end  of  the  body  will  all  originate 
near  the  wine  pV.-ie  at  the  initial  plane. 

The  details  on  the  streamline  calculation  procedures  used  on 
the  nosetip  and  the  frustum  are  provided  in  Sections  3.3.1  and  3.3.2, 
respectively,  below. 

3.3.1  SPHERICAL  NOSETIP  STREAMLINE  CALCULATION 

The  relationship  between  the  3IS  (r,0,z)  cylindrical  coordinate 
system  and  the  2IT  (R,  O',  V)  wind-fixed  spherical  coordinate  system  is 
depicted  In  Figure  3.1.  The  locations  of  the  required  inviscid  surface 
streamlines  at  the  3IS  initial  data  plane  are  assumed  to  be  known,  and 
are  defined  by  specifying  the  appropriate  value  of  9  for  each  streamline. 
The  track  of  each  of  these  streamlines  on  the  nosetip  is  readily  defined 
by  0'  =  constant  on  the  sphere.  The  problem  thus  reduces  to  determining 
the  relationship  between  9  (3IS  coordinates)  and  o'  (2IT  coordinates)  at 
the  initial  data  plane. 
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21 i  and  31S  Coordinate  Systems 


where  G^  Is  the  body  angle  at  the  Initial  data  plane,  is  the  nose 
radius,  and  z^  is  the  distance  of  the  initial  data  plane  from  the  virtual 
apex. 

Knowing  the  appropriate  value  of  o'  for  each  desired  streamline 
the  inviscid  nosetip  flow  field  parameters  may  then  be  defined  at  the 
desired  axial  stations  (in  3IS  coordinates),  requiring  only  transforming 
the  2IT  velocity  components  to  the  appropriate  3IS  velocity  components. 
Given  the  U, V  velocity  components  of  the  2IT  solution  (defined  in 
References  2  and  3),  the  3IS  velocity  components  (u,v,w)  may  be  written 
at  a  given  station  (z)  as 


u  =  -V  [s  1  n<j>  cose  sina  +  coscj)  (sine  sine'  +  cose  cose'  sina)] 

+  U[-cos((>  cosO  sina  +  sincf>  (sine  sinG*  +  cose  cose'  cosa)] 


(3.2 


v  =  -V[-sin<j>  sinO  sina  +  cos<i>  (cose  sine'  -  sine  cose'  cosa)] 
+  U [cos4>  sine  s1na  +  s1nt)>  (cose  sine'  -  sine  cose'  cosa)] 


(3.3; 


w  -  -V  (sinc|)  cosa  -  cost))  cose'  sina) 
-U  (sini|i  cose'  sina  +  cos<b  cosa) 


(3.5' 


where 


cos  <t>  E  cosa  cosv  -  sina  sinv  cosG 


cos  v  =  1 


Z  -  Zp 

Rn 


tan 


_ _  sinij)  sine* _ 

sTrvfi  cosOT  cosa  -  cosep  sina 


with  z0  being  the  intersection  of  the  3IS  axis  with  the  spherical  nose. 


The  nosetip  flow  field  data  obtained  using  this  procedure  is 


structured  in  the  format  described  in  Section  3.4.2  for  use  by  the  3DMEIT 
code. 

3.3.2  AFTERBODY  STREAMLINE  CALCULATION 

On  the  afterbody,  as  noted  earlier,  a  procedure  has  been  de¬ 
veloped  in  this  effort  to  allow  the  explicit  calculation  of  inviscid 
surface  streamlines  within  the  3IS  code.  This  procedure  takes  advantage 
of  the  fact  that  functions  which  are  conserved  along  streamlines  can  be 
readily  computed  in  an  inviscid  flow  field  code  using  convective  finite 
difference  schemes.  The  calculation  of  entropy  in  the  3IS  code  is  an 
example  of  such  a  procedure.  (Since  inviscid  flows  are  isentropic,  the 
value  of  entropy  does  not  change  along  a  streamline  in  Inviscid  flows.) 

The  attributes  of  such  "convective"  functions  can  readily  be 
used  to  simplify  the  calculation  of  Inviscid  surface  streamlines,  In 
the  3IS  code,  each  surface  streamline  is  identified  by  a  unique  value  of 
a  function  f,  which  is  chosen  to  represent  the  value  of  9  in  the  Initial 
data  plane  through  which  that  streamline  passes.  The  "conservation  of  f" 
along  streamlines  can  then  be  expressed  mathematically  as 

•  V f  =  0  (3.5 

which  can  be  expressed  In  terms  of  the  3IS  variables  at  the  body  surface 
as 


9f  _ 
9Z 


(3.6 


where  Y  is  the  transformed  circumferential  coordinate  and 


(For  more  detail  on  the  3IS  equations  and  transformations,  see  References 
2  and  3.) 

Equation  (3.6)  is  then  solved  using  second-order  accurate  differences  in 
a  predictor-corrector  scheme,  simultaneously  with  the  other  flow  equations 
being  solved  in  3IS,  to  define  the  appropriate  value  of  f  at  all  grid 
points  located  on  the  body. 

Given  the  computed  f  distribution  over  the  entire  afterbody 
surface,  the  track  of  any  given  streamline  (f  =  constant)  can  rapidly 
be  determined  by  Interpolation,  along  with  the  rest  of  the  flow  field  in¬ 
formation  to  be  passed  to  3DMEIT  by  the  buffer  code. 

3.4  USER  INSTRUCTIONS 

3.4.1  INPUT  AND  OUTPUT 

As  input,  the  buffer  code  between  the  3IS  and  3DMEIT  procedures 
requires  only  two  parameters  to  be  Input:  the  nose  radius  (RN)  and  the 
body  angle  at  the  nosetip  tangency  point  (THB1),  in  degrees.  These 
parameters  are  input  in  NAMELIST  format.  All  other  required  inputs  are 
provided  either  through  the  2IT  output  file  (which  must  be  assigned  to 
logical  unit  23)  and  the  3IS  output  file  (which  must  be  assigned  to 
logical  unit  21).  The  formats  of  the  2IT  and  3IS  output  files  are 
described  in  Reference  3. 

The  output  of  the  buffer  code  is  written  to  a  formatted  file 
assigned  to  logical  unit  22.  The  format  of  this  file,  which  will  be  read 
as  input  to  the  3DMEIT  code,  is  described  in  the  following  section. 


20. 


3.4.2 


INVISCID  FLOW  FIELD  FILE  DESCRIPTION 


tains: 


The  first  record  on  the  output  file  of  the  buffer  code  con- 


Mach  number 


a  Angle  of  attack 

MMAX  Number  of  streamlines 

NMAX  Number  of  points  through  the  shock  layer. 

The  format  of  this  record  is  2E13 .6,213 ,48X . 

The  second  record  contains  the  axial  location  on  the  body,  Z, 
in  the  format  E13.6.67X. 

The  third  record  contains  values  for 
0  Circumferential  angle 

rb  Body  radius 

rbz  3rb/3z 

rb0  3rb/30 

oj  tan_1(v/w),  where  v  and  w  are  the  circumferential 

and  axial  velocity  components,  respectively 

SF  Stream  function,  used  to  identify  individual 

streamlines 

for  each  of  the  MMAX  streamlines  at  the  axial  station  defined  in  the 
preceding  record.  These  data  are  written  in  the  format  6E13.6.2X. 

The  next  NMAX  records  contain  values  for 
p  Pressure 

T  Temperature 

u  Radial  velocity  component 

v  Circumferential  velocity  component 


2 


w  Axial  velocity  component 

y  Distance  from  the  body  surface,  measured 

radially 

for  each  streamline,  written  in  the  format  6E13.6.2X.  Each  of  these 
records  corresponds  to  a  different  value  of  N,  where  N  =  1  corresponds 
to  the  body  surfnce  and  N  =  NMAX  corresponds  to  the  bow  shock. 

This  sequence  of  recorJs,  starting  with  the  second  record,  is 
repeated  for  each  axial  station  at  which  the  3IS  solution  produced  complete 
field  output. 
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SECTION  4 


BOUNDARY  LAYER  SOLUTION 

4.1  INTRODUCTION 

The  3DME IT  code  is  a  solution  procedure  for  the  boundary  layer 
integral  momentum  and  energy  equations  in  a  streamline-body  normal  ’• 

coordinate  system  over  a  general  three-dimensional  body.  The  flow  may  be 
laminar,  transitional  or  turbulent,  and  the  procedure  accounts  for 
compressibility,  real  gas  effects  and  surface  ablation.  The  usual  inputs 
to  the  code  are  the  geometry  of  the  body  and  the  streamlines,  a 
specification  of  inviscid  flowfield  data  near  the  body  surface  and 
selection  of  a  wall  boundary  condition. 

Closure  of  the  set  of  integral  equations  is  accomplished  by 
specifying  the  local  shape  factor,  the  recovery  factor,  the  Stanton  number 
and  the  friction  coefficient  as  functions  of  the  momentum  and  energy 
thickness  Reynolds  numbers.  These  basic  "laws"  are  modified  to  account 
for  the  effects  of  surface  roughness,  transpiration,  acceleration,  and 
compressibility  by  influence  coefficients,  which  are  multiplicative 
factors  on  the  local  Stanton  number  and  friction  laws.  The  numerical 
solution  procedure  is  an  implicit  finite  difference  scheme. 

In  Section  4.2,1  the  basic  equations  are  presented.  This  is  followed 
by  the  formulation  of  local  shape  factor,  recovery  factor,  Stanton  number, 
friction  coefficient  and  influence  coefficients  in  Section  4.2.2,  The 
transition  criteria  and  surface  roughness  modeling  are  presented  in 
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Sections  4.2.3  and  4.2.4,  respectively.  Section  4.2.5  describes  the  models 
used  to  approximate  the  ablation  effects  while  Sections  4.2.6  and  4.2.7 
detail  the  entropy  swallowing  and  induced  pressure  models.  The  solution 
procedure  is  described  in  Section  4.2.8 

4.2  3DMEIT  METHODOLOGY 

4.2.1  Basic  Equations 

The  two  equations  that  are  basic  t.o  3DME1T  are  the  integral 
momentum  equation, 


ae  _ H _ 8£  1_  _ 1 _ (  2v  +  £f  .  pwvwuiw 

’S  "  “eue2  "  “  "  peue2  “  *  6  2  'peVe 


(4.1) 


and  the  integral  energy  equation. 
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where  s  is  distance  along  a  streamline  and  h0  is  the  body  surface 
metric.  These  forms  of  the  boundary  layer  equations  result  from 
neglecting  the  crossflow  velocity;  3DMEIT  is  therefore  most  appropriate 
for  regions  where  the  crossflow  is  "small."  However,  a  simplified  form  of 
the  crossflow  momentum  equation  can  be  added  to  3DMEIT  without  changing 
the  parabolic  nature  of  the  system  of  equations.  This  addition  is  a 
relatively  simple  change  to  the  code  and  should  be  considered  if 
predictions  at  the  next  level  of  detail  in  the  flow  are  necessary. 

The  other  variables  in  these  equations  ere  the  momentum  and  energy 


thicknesses,  which  are  respectively: 


J  peue\  ue  / 


(4.3) 


no 

=  /(5 


_  hte  -  ht 
e  ^te  "  hw 


(4.4) 


and  the  boundary  layer  shape  factor,  H,  is  defined  as: 


where  6*,  the  displacement  thickness,  is  given  by: 


(4.5) 


•/(■•*) 


(4.6) 


The  total  enthalpy  at  the  boundary  layer  edge  is  defined  by: 


hte  "  he  +  r 


(4.7) 


while  the  recovery  enthalpy  is  given  by: 


hy.  =  h  +  F  -5-- 

re  £ 


(4.8) 


where  F  is  the  recovery  factor. 

The  heat  transfer  rate  and  skin  friction  are  related  to  the  Stanton 
number  and  friction  coefficient,  respectively,  by: 


Tw  =  peue  T 


(4,9) 
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Closure  of  this  system  of  equations  is  described  in  the  following 
subsection. 

4.2.2  SHAPE  FACTOR,  RECOVERY  FACTOR,  STANTON  NUMBER  AND  FRICTION 
COEFFICIENT 

The  shape  and  recovery  factors  are  evaluated  in  3UMEIT  by  the 
following  relation  for  laminar  flow: 

H0  -  3.029  I?-  -  0.614  (4. 

*  'e 

F£  =  Pr1/2  (4. 


and  for  turbulent  flow: 

Ht  =  Href  [1  +  D(6  -  log10Ree)] 


(4. 


where 


Href  =  A  e 


BMl 


A  = 


Ao  +  A1 


0.0064  +  0.0231  (Tw/Tte) 


IK 


A0  .  -0.02  +  1.20  (Tw/Tte) 
A1  =  (0.0372  +  0.0322 


(4. 

(4. 

(4. 

(4. 
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log10(ks/3) 


0.35  (VTter0,5 
3.817 


for  ks  <  3 
for  >  3 


for  M  <  1.5 
for  M  >  1.5 


for  M  <  1.5 
for  M  >  1.5 


(4.18) 


(4.19) 


(4.20) 


Mm  M 


for  k  <  3 


0.05  +  0.04  log10(ks/3)  for  k  s  >  3 


(4.21) 


(4.22) 


Ft  =  PrJ 


(4.23) 


where  kg  is  the  local  roughness  Reynolds  number. 

The  friction  coefficient  and  the  Stanton  number  are  determined  by 
the  basic  wall  shear  and  heat  flux  laws  respecti vely.  The  basic  form  of 
these  laws  comes  from  theory  and  data  for  incompressible  flow  along  a 
smooth,  isothermal,  impervious,  flat  plate.  The  friction  coefficient  and 
Stanton  number  for  laminar  flow  are: 


'f.i.o  0.245 

~2  =  “ReT- 


(4.24) 


V 


(4.25) 


M,o  ^ 


and  for  turbulent  flow: 


Cr  j.  .  /->  r  0.010742  Re.  .  coco 

T,t,o  0.24b  A  e  n _  \-1.5262 

- re—  +  T0r+_R''eT'  (,og10  Ree} 

e  e 


(4.26) 


0.22  + 


^,t,o  pr4/3Re^  (ioo  +  Re^) 


;rr  °09io  V6 


(4.27) 


where 


a  -  0.0993  e(0,0648  Pr^/(1  +  9.6  Pr) 
b  .  1.954  e(0‘273  Pr)/(1  +  0.71  Pr) 


Pr  >  0.6 


(4.28) 


a  -  0.1256  e(0,2435  Pr)/(1  +  14.2  Pr) 
b  -  2.217  e(°*6313  p r)/(!  +  1.677  pr) 


Pr  <  0.6 


(4.29) 


The  Stanton  number  and  friction  coefficient  laws  given  above  are 
modified  by  corresponding  influence  coefficients  to  account  for  various 
boundary  layer  effects: 


z 

c  .  c  n  i 

x,y  x,y,o  x,y,z 


x  =  h,f 
y  «  M 


(4.30) 


The  influence  coefficients  are  shown  by  I  where  the 

*  »y  •  * 

subscripts  x  and  y  indicate  whether  the  influence  coefficient  pertains  to 
heat  or  momentum  transfer  (x  =  h  or  f)  and  laminar  or  turbulent  flow 
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(y  ■  1  or  t),  respectively.  The  subscript  z  Indicates  the  type  of 
phenomenon  for  which  the  basic  laws  are  being  modified. 

Five  phenomena  are  considered  by  3DMEIT.  These  phenomena  and  their 
corresponding  z-subscrlpts  are  given  in  Table  4.1.  The  influence 
coefficients  corresponding  to  each  of  these  effects  are  given  below. 
Acceleration 

In  laminar  flow: 


(1 + 

a  >  o 

1.0, 

a  <  o 

(1 + 

4«)  . 

a  >  o 

1.0, 

a  <  o 

(4.31 


(4.32 


where 


B  _  25  dUe 


(4.33 


■  /  p<=" 


au.r  ds 
e  e 
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Table  4.1.  Notation  for  Influence  Coefficients 


Phenomena 


z-Subscript  of  Influence 
Coefficient 


Acceleration 
Transpiration 
Boundary  layer  properties 
Roughness 

Transition  proximity 


In  turbulent  flow,  only  the  first  terms  of  the  appropriate  basic 


turbulent  laws  are  modified  by  the  above  influence  coefficients. 
Transpiration 

Blowing  effects  on  the  wall  shear  are  modeled  by  adaptations  of 
film  theory  (Reference  14). 


M  i  +  2xf 


!f ,y,B'  "  2xf  *  v  "  A»t 


(4.35) 


where: 


B'  = 


<n/0«u_CL 
e  e  m 


(4.36) 


and 


R  -  Ch/Cf/2  -  Reynolds  Analogy  Factor 

R'  -  Mass  to  heat  transfer  coefficient  ratio  =  Cm/C^ 

X..  •  Blowing  reduction  parameter. 

■  *y 

R  is  a  dependent  variable  which  is  evaluated  during  the  solution 
process.  Both  R'  and  xf  are  input;  3DMEIT,  however,  does  provide 
built-in  default  values  for  xf  .  These  default  values  are: 

xf,A"  0,5 

(4.37) 

^f ,t  “  0*35 

The  variables  B'  and  R'  are,  in  effect,  definitions  for  the  mass 
transfer  coefficient,  which,  in  3DMEIT,are  obtained  from  the  analogy  be¬ 
tween  heat  and  mass  transfer.  Therefore,  Cm  is  the  mass  transfer  analog 
of  the  Stanton  number  and  m  Peue^m’ 
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Ordinarily,  R1  s  1  and  the  default  values  of  A  are  used  because 
these  are  the  theoretical  values  for  unity  Prandtl  and  Lewis  numbers  in 
laminar  flow.  Also,  the  default  values  for  turbulent  flow  have  been  chosen 
to  fit  available  data.  Other  input  values  of  R1  and  A  can  be  developed 
from  finite  difference  calculations  for  chemically  reacting  boundary  layers 
or  by  fits  to  mass  transfer  data.  (A  detailed  description  of  these  pro¬ 
cedures  is  given  in  "An  Evaluation  of  a  Transfer  Coefficient  Approach 
for  Unequal  Diffusion  Coefficients"  by  E.  P.  Bartlett  and  R.  G.  Grose, 
Aerotherm  Report  69-50,  June  30,  1969.) 

There  are  two  blowing  correlation  options  for  the  heat  transfer: 

(1)  Aerotherm  option  and  (2)  GE  option.  The  Aerotherm  option  utilizes  the 
expression  (Reference  15): 


1 H ,y , B 1 


y  -  *.t 


(4.38) 


The  GE  option  employs  the  Costello  correlation  (Reference  16)  for  laminar 
flow: 


'h.n.B'  -tr-  <»  ‘O^BT1 


(4.39) 


and  another  correlation  (Reference  17)  for  turbulent  flow 


1 h, t , B ' 


(1  +  B') 


0. 2  ( ui— 1 ) 


(4.40  ) 


where 


+  1/8  M, 


(4.41) 


Boundary  Layer  Properties  . 

Boundary  layer  properties  of  density,  viscosity,  and  Prandtl  number 
are  evaluated  at  a  reference  enthalpy  h‘  as  follows  for  laminar  flow: 


h'  -  0.42  hfi  +  0.19  hr  +  0.58  hw 


(4  42 ) 


and  the  property  influence  coefficients  are: 


'f.t.p  ■  1 


(4.43) 


M,p 


^p'\  / u 1  \  /hte_  “  hw\ 


(4.44) 


For  turbulent  flows: 


h'  -  0.353  he  +  0.19  hr  +  0.45  hw  +  fg  f(he,hr,hw)  (4.45) 


where 


t(hg,hr,hw)  is  the  greater  of  zero  or  (-0.609  h£  +  0.332  hr  +  0.277  h  ) 


(4.45) 


(kg  -  3)2 
400  +  (k*  -  3)‘ 


for  ks  <  3 


for  k$  >  3 


(4.47) 
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where  k*  is  the  local  roughness  Reynolds  number. 

The  influence  coefficients  are: 

\ 

/  ,  v  /  ,  ,0.25 

Ix*t.P  ’  X  "  f>h 

Surface  Roughness 

(4.4B 

In  laminar  flow: 

1  h, JL . r  *  *f,£tr  *  ■ 

(4.49) 

The  turbulent  influence  coefficient  for  the  wall 

shear  due  to  roughness  is: 

a(log10  Re,)-e 

tA  t:n  \ 

0.010742  (log10Ree) 


where 


a  =  b  +  c(k*)d 


e  -  f  +  g( ) h 


The  values  of  b,  c,  c,  f,  g  and  h  are  dependent  on  the  local 
roughness  Reynolds  number  (kj  a  uTks/vw)  and  are  listed  in 
Table  4.2.  Ree  the  local  momentum  thickness  Reynolds  number  for  the 

rough  wall,  uT  is  the  local  shear  velocity  (uT  A  s/r^o^),  ks  is  the 

equivalent  sand-grain  roughness  height  and  v  is  the  wall  kinematic 

w 

viscosity. 

The  turbulent  heat  transfer  influence  coefficient  is: 


R  -  fl  [1  +  92<93  -  f2)][1  -  f3]  +  f3  “  91 


h 


(4.54) 


*e 


Table  4.2, 


Coefficients  for  the  Turbulent  Wall  Shear  Roughness 
Influence  Coefficient 


[log10( 10,000  ks/e)2 J  for  ks/e  >  10"4 
0  for  ks/e  £  10"4 

gj  *  0.005  {1  +  0.35  sin  [vr/2  (1  -  Pr)]f 
g2  =  0.142  (Pr  -  0.71) 

g3  .  (0.121  +  0.053  Pr)-1 

The  values  f 2  and  f^  are  functions  of  k*  and  are  given  in 
Table  4.3. 


(4.55  ) 

(4.56) 

(4.57) 

(4.58) 
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Table  4.3.  Functions  f2  and  f 3  for  the  Turbulent  Heat 
Transfer  Roughness  Influence  Coefficient 


Range  of  k* 

f2 

f3 

3  -  76 

1.8808 

-0.32  [  log1Q(kg  -  1.8808)] 

76  -  1500 

1°gio<ks) 

0 

1500  -  104 

3.176 

0.32  [log1Q(kg  -  3.176)] 

The  function  f^  is  the  relationship  between  heat  transfer  and 
roughness  shear  augmentation.  The  modified  correlation  is: 


0.89  -  0.46<k*  -  3) /{k^  +  47) 
0.89 


for  >  3 
for  kg  <  3 


(4.59  ) 


Transition  Proximity 

The  transition  proximity  influence  coefficients  apply  only  to  the 
laminar  region  upstream  of  the  transition  point.  The  present  influence 
coefficients  are  based  on  PANT  and  ART  wind  tunnel  heat  transfer  data 
(References  18  and  19).  It  is  assumed  that  the  influence  coefficient  for 
momentum  is  the  same  as  for  heat  transfer,  i.e., 


Jf, Jl.tr  “  !h, Jl.tr 


(4.60) 


and 


Ih  ^  tr  =  the  greater  of  1.0  or 


7  ,,  ,-0.366  °*0005  stk  ,nnn 

7 . 6  ( s  t  k )  e  s+l.  <  1000 


*tk 


(4.61) 
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where 


stk 


str/etr^h,il,tr^ 


1.3 


(4.62) 


Intermittency 

The  above  formulations  for  H,  F,  Ch,  C are  for  entirely 
laminar  or  for  entirely  turbulent  flows.  To  evaluate  these  four 
parameters  for  transitional  flow  the  following  relation  is  used: 

P  «  (1  -  f)  P*  +  fPt 


where  P  is  one  of  the  four  parameters  above  and  f  is  the  transitional 
intermittency  factor. 

The  transitional  intermittency  employed  in  3DMEIT  is  based  on  the 
work  of  Persh  (Reference  20),  and  an  interpretation  by  Dahm  (Reference  21): 


f  -  1 


Ree  (Cf,t  '  Ct'A 


(4.63) 


where 


a 


Ree,tr  (Cf,t  "  Cf,Ar 


(4.64) 


and  the  subscript  tr  refers  to  conditions  at  the  transition  point. 

f  is  set  to  zero  in  laminar  flow,  unity  in  turbulent  flow,  and 
varies  between  0  and  1  in  transitional  flow. 

This  completes  the  formulation  of  the  parameters  required  in  the 
solutions  of  Equations  (4.1)  and  (4.2). 
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The  3DMEIT  code  has  several  built-in  options  to  specify  the 


critical  flow/material  surface  parameters  that  establish  transition  onset 
and  location.  In  addition,  all  laminar  or  all  turbulent  (from  stagnation 
point)  calculations  may  be  made. 

The  analytic  transition  criteria  for  rough  wall  calculations 
available  are  the  Anderson  criterion  for  nose  transition  and  LORN  criterion 
for  cone  transition  (although  use  of  the  LORN  criterion  is  not  recommended) 
These  are: 

Anderson  Criterion 


255  onset 
215  location 


(4.65 


l 

where 


\l>  =  0.1B'  +  (1  +  0.25  B')  p6/pw 


(4.66 


LORN  Criterion 

0.134  M 

Re  =  275  e 
e 


(4.67 


where  Mg  is  the  boundary  layer  edge  Mach  number. 

4.2.4  SURFACE  ROUGHNESS  MODELING 

Two  types  of  surface  roughness  are  modeled 

•  Intrinsic  roughness,  k- 

•  Turbulent  or  scallop  roughness,  kt> 

Intrinsic  roughness  (k^)  is  that  associated  with  the  basic  material 
granularity  and  is  input  as  a  constant  for  each  material  or  as  a  function 


of  the  axial  length.  The  intrinsic  roughness  is  used  In  the  transition 
criteria  and  transition  proximity  influence  coefficient. 

The  turbulent  roughness  (k^)  is  the  effective  sand  grain 
roughness  that  results  from  turbulent  ablation  and  is  used  in  the 
turbulent  roughness  influence  coefficients.  The  roughness  height,  kt, 
is  specified  in  one  of  three  ways.  A  uniform  value  of  kt  for  all 
turbulent  regions  may  be  input  by  the  user;  a  distribution  of  k^.  as  a 
function  of  axial  length  may  be  input;  or  a  value  may  be  obtained  by  using 
the  scallop  dimension  correlation.  From  the  correlation,  the  effective 
turbulent  region  scallop  depth  is  computed  as  follows: 

kt  =  Kipe~°‘77  (4.68) 

where 

kt  -  the  effective  sand  grain  roughness  height 

Kj  -  a  material  dependent  property  determined  from  experimental 
data  and  input  by  the  user. 

4.2.o  SURFACE  TEMPERATURE  AND  ABLATION  MODELING 

The  technique  used  in  3DMEIT  is  an  approximation  because  it  assumes 
steady-state  conditions;  i.e.,  no  transient  conduction  in  the  solid. 

Under  these  special  conditions  of  quasi-equilibrium,  the  energy  balance  is 
(excluding  mechanical  erosion): 


qNET  ”  qc  “  qR 


mxhlf  +  m,h  +  m„h„ 
T  w  c  c  g  g 


(4.69  ) 


where 


m-r 


m. 


(4.70) 


■i 


I 

"i 
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Since  the  net  heat  conducted  into  the  solid  is  approximately  the  energy 
lost  due  to  mass  removal,  1 . e. : 


qNET  =  \\  *  ligh,  (4.71) 

The  surface  energy  balance  can  be  further  simplified  to  give: 

De"eCh(hr  -  hw>  -  "Th»  -  '  0  <«•»> 


This  equation  Is  solved  iteratively  for  and  hw  by  varying  the  wall 
temperature  until  the  recovery  enthalpy  matches,  within  a  specified  limit, 
the  local  value  derived  from  mass  balance  considerations.  The  above 
equation  is  conveniently  rewritten  for  this  purpose  as: 


ah„  =  1.0  - 
r 


Ca1w  +  mThjw  + 


^eue^h 


V/hr 


(4.73) 


4. 2. 5.1  Hall  Gas  Enthalpy 

The  wall  gas  enthalpy  for  graphitic  materials  is  expressed  in  the 
functional  form: 

hw  -  f  (B‘)  (4.74) 

where  B'  is  the  blowing  parameter: 


B' 


”T/'e“eCh 


(4.75) 


and  Ch  is  the  Stanton  number,  Pconvective/peue^hr  ~  hw)»  with  blow1n9* 

f(B')  is  represented  by  the  following  polynomial  relationships  for 
the  range  of  indicated  pressures: 


X. 


A  -  1.4216  +  0.0293  log1Q  pg 
C  -  0.55993  +  0.00348  logjo  Pe 

(4.76) 

D  -  -8.1997  -  0.2173  log1Q  pg 
E  =  0.9024  -  0.062  log1Q  pg 


B'  <  0.19:  h*  -  [-  E  +  \/  E  2-  4C  (0.0 361 A  +  0.19D  +1)  ] 

hwD  -  -172.46  +  76.63  log1()Pe  +  20.41  (log10Pe)Z  (4.78) 
h»  ■  [<hwJ  -  hw  >  B'  +  Vs  hw  -  °-m  K  ]  '  °-015(4.79) 


B1  >  0.19: 


2500 

2C 


-  E  +  n/e2  -  4C  (AB|2  +  DB*  +  1) 


(4.80) 


These  relationships  also  satisfactorily  approximate  the  wall  gas 
enthalpies  for  carbon  phenolic  (Reference  15) . 

A  discussion  of  the  relevant  expressions  used  in  the  three  ablative 
material  options  is  presented  in  the  following  subsections: 

4 . 2 . 5 . 2  Graphite  Ablation  Option 

For  graphite  (carbon/carbon,  ATJ,  etc.)  ablation,  the  phenolic 


resin  gas  flowing  rate  is  zero;  i.e.,  mg  -  0.  The  ablation  rate 
expression  employed  is  of  the  form: 


B1  =  B*  =  exp 

V 


5 


2 


n 

bnT 


(4.81) 
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where 


4 

T*  '  Tw  -  6000  -  X  \  ’°9l0  Pe 

n«l 

and  the  numerical  values  for  a„  and  b  are  listed  as  follows: 

n  n 


Ablation  Parameter  Correlation  Constants 


n 

an 

bn 

0 

- 

-1.3191 

1 

587.496 

6.5704  (10“4) 

2 

32.826 

3.6955  (10~7) 

3 

9.643 

6.7268  (10"1Z 

4 

-2.965 

-2.7179  (10-14 

5 

3.3450  (10“17 

T*  is  employed  as  a  collapsing  function  which  reduces  the  data  to  one 
curve.  The  relationship  between  T  and  B^  is: 

T*  <  -4447:  *  0.175 

..4447  <  T*  <  -939:  8^.  -  0.175  +  4.2759  ( 106)  (T*  +  4447) 

T*  >  -939:  =  exp 


5 

2 

Ln-0 


b_  T 


(4.82) 


(4.83) 


41. 


4 . 2 . 5 . 3  Carbon  Phenolic  Ablation  Option 

In  the  carbon  phenolic  option,  the  total  mass  loss  rate  from 
Equation  (4.70)  becomes: 


where  the  char  mass  loss  rate  parameter  is  obtained  from  Equations 
(4.81)  and(4.82)  and  the  ratio  /mc  (correlated  as  a  function  of  the 
zero-blowing  film  coefficient  Ch0  from  results  reported  in  Reference  22  )  is 


1 


0*22  Ch  ^  j  q 

0.22  -  0.01  Ch  Ch  >  1>0 


(4.85) 


The  wall  gas  enthalpy  for  carbon  phenolic  ablation  Is  satisfactorily 
approximated  (see  Reference  23 )  using  Equation  (4.74)  in  the  form 


hw  ■  f<e;> 


(4.86) 


This  approximation  is  valid  only  for  mg/mc  <  10  which  encompasses 
most  reentry  applications. 

The  procedure  to  compute  the  blowing  rate  and  wall  temperature  is 
identical  to  that  described  in  the  "Graphite  Ablation  Option"  section. 

4 . 2 . 5 . 4  Teflon  Ablation  Option 

For  Teflon  ablation,  the  empirical  "heat  of  ablation"  concept  is 
employed  to  obtain  the  mass  loss  rate.  The  expression  employed  in 


Til 


'll 


i 


0  ’ 


■'  1 
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subroutine  ABLATE  was  obtained  from  Reference  24  and  is  utilized  in  the 
form: 

ch 

L  0 

mT  "  074 b 94  -  0.0027  pe  (4.37) 


with  the  ablation  temperature  of  Teflon  set  at  1,800°R.  Any  of  the 
three  material  options  described  above  may  be  used  at  one  time. 

4.2.6  EDGE  ENTROPY  MODELS 

The  3DMEIT  code  has  four  different  models  for  determining  the 
entropy  at  the  edge  of  the  boundary  layer.  The  first  two  options  are  very 
simple  and  do  not  account  for  entrainment  into  the  boundary  layer. 

Option  1  uses  a  constant  entropy  equal  to  that  behind  a  normal  shock  and 
Option  2  uses  the  entropy  calculated  from  the  invlscid  pressure  and 
temperaure  at  the  wall.  Options  3  and  4  account  for  swallowing  and  use 
the  invlscid  flowfield  data  to  calculate  the  edge  entropy.  Option  3  Is 
the  streamtube  entrainment  method  used  in  ASCC  and  Is  shown  schematically 
in  Figure  4.1, 

Balancing  the  mass  entering  the  streamtube  with  the  mass  flux  in 
the  boundary  layer  gives: 


/ 


p»Vdy  "  2hB  I  pU  dy  _  2 


/' 

HP 


hB(pV!wdS 


(4.88) 


o'  o*' 

where  hg  is  the  spreading  metric.  In  axi symmetric  flow,  hg  =  r(s),  which  is 
the  local  body  radius.  Note  that  y  is  the  body  normal  coordinate  in  the  right 
hand  side  of  this  equation.  Then,  using  MEIT's  variables,  this  becomes: 


p  u  y 

00  co*' 


f  hB(pV>wdS 

r 


2hp  Fue  Re0  -  2 


(4.89) 


Figure  4.1.  Schematic  of  Boundary  Layer  Entrainment 


Equation  (4.89)  Is  solved  for  y  which  Is  then  used  with  the 
Inviscid  shock  shape  and  entropy  to  determine  the  entropy  of  the 
streamline  at  the  edge  to  the  boundary  layer.  This  model  works  well  on 
ax  i  symmetric  configurations  at  zero  angle  of  attack.  However,  in 
three-dimensional  problems,  the  tracking  of  the  streamline  at  the  edge  of 
the  boundary  layer  is  difficult  because  the  streamlines  curve  around  the 
body.  Therefore,  Option  4  was  developed  to  utilize  a  mass  balance  between 
the  inviscid  flow,  as  defined  by  a  finite  difference  calculation,  and  the 
boundary  layer  flow  (Figure  4.2).  Balancing  the  mass  fluxes  at  a  local 
station  gives  for  this  case: 


/  <»u>1nv  dy  '  /  t>u)  BL  dy  -  / 


BL  dy  -  /  ha(Pv)wds  (4.90) 

_  J  J 

0  0  0 

where  6  is  used  as  the  integration  limit  for  both  the  inviscid  and  boundary 
profiles  because  the  "edge"  of  the  boundary  layer  must  have  the  same  p  and  u 
as  that  which  exists  in  the  inviscid  flow  at  that  6  and  the  boundary  layer  mass 
flux  must  Include  that  which  is  entrained  from  the  entropy  layer.  This  equation 
effectively  defines  5,  which  is  otherwise  unknown  and  not  needed  for  an 


integral  boundary  layer  calculation 

f  (»u))nv  <<y  -  »,R«eF  -  f  V»v)wd5 


(4.9l) 
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The  right  hand  side  of  this  equation  Is  known  from  the  boundary  layer 
solution.  The  solution  for  6  is  obtained  by  constructing  the  inviscid 
mass  flux  as  a  function  of  the  distance  normal  to  the  surface  and 
interpolating  In  these  inviscid  data. 

4.2.7  INDUCED  PRESSURE 

The  effect  of  the  boundary  displacement  on  the  pressure  is  modeled 
by  a  correction  to  the  inviscid  values 


P 


p, 

Hinv  ae  dx 


(4.92) 


where  da*/dx  is  the  derivative  of  the  displacement  thickness  with  respect 
to  the  axial  direction.  The  term  ap/ab  models  the  effect  of  the  angle  on 
the  wall  pressure.  This  term  is  evaluated  using  Linnell's  formula 
(Reference  25)  which  may  be  written  as: 

jj-  -  1  +  0.5  Y  (y  +  1)  K2  (4.931 

*00 

where  K  is  the  hypersonic  similarity  parameter,  K  -  sin  e. 

The  term  ap/ae  is  determined  by  calculating  p  at  the  body  angle  and 
at  the  effective  angle: 

e'  =  e  +  tan-*  (as*/as)  (4.94, 

When  e  is  less  than  zero,  limiting  values  of  e  «  1°  and  e*  »  2°  are  used 
in  calculating  ap/ab. 

4.2.8  SOLUTIONS  OF  THE  BOUNDARY  LAYER  INTEGRAL  EQUATIONS 


As  mentioned  above,  the  required  input  to  3DMEIT  are  surface  shape, 
boundary  layer  edge  conditions,  and  wall  conditions.  These  quantities  are 
input  in  terms  of  body  points  and  inviscid  solution  points.  A  finer  grid 


in  terms  of  integration  points,  which  include  all  the  inviscid  points,  is 
generated  by  the  program  to  ensure  adequate  solution  accuracy  of  the 
integral  equations.  The  boundary  layer  edge  conditions,  gas  properties 
and  wall  conditions  at  each  integration  point  arc  obtained  by 
interpolation  from  the  input.  The  solution  procedures  of  the  boundary 
layer  integral  Equations  (4.1)  and  (4.2)  consist  of: 

•  Start-up  solutions  at  the  first  integration  points 

•  Finite  difference  numerical  solutions  for  the  rest  of  the 
Integration  points. 

The  starting  values  of  <4  and  e  are  defined  in  the  code  by: 


el  " 


0.245  vj  (1  '+  R 

w 

du 

<3  +  H)  is" 

1 

77  C 


f.Jt.2 


(4.95) 


0.22  v-  (1  +  Bt)  Z 

1  77  C 


2  Pr4^3  ^ 


M,z 


(4.96) 


The  solutions  at  the  second  and  third  integration  points  are 
related  to  the  first  integration  point  by: 


"'ll1*  a*2) 

(4.97) 

*  ■  <1  (1  * 

(4.98) 
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where 


13  +  H1  a  .  1  °‘6S9  <3  +  Hl>  a  (>  -  1)  (H.,  -  0.614)  a  (y  -  1) 

4 -  y'? - - - - - - - 


8 


(4.99) 


1  [  3a  +1_  ?  X  0.659  a  (y  -  1|  2  ^ 

6  l  y  3  y 


-  F.)(]  +  B*  )  a  (y 


(4.100) 


*  “  S/Rref 


(4.101) 


i  -  ( P3/Pi ) 
~2 


(4.102) 


In  the  above  formulation,  R  ^  is  an  arbitrary  constant  radius  and 
y  is  the  isentropic  exponent.  The  subscripts  1  and  3  denote  the  first  and  third 
integration  point  condition,  respectively  (e.g.,  Hx  is  the  shape  factor  at  point  1). 

The  following  implicit  finite  difference  scheme  is  used  for  the 
rest  of  the  integration  points: 


Fx,I  ■  Fx,I-l  +  °-5  (fi.i-i  *  -  *!_!>  X  -  f.h 


(4.103) 


where: 


Ff  *  r  peuete 


[ see  Equation  (4.1)] 


(4.104) 


Fh  "  r  peue^hte 


[ see  Equation  (4.2)1 


(4.105) 


>  *  Us” 


(4.106) 
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(4.107) 


F 


h 


dFh 

nr 


I  Is  the  Integration  point  Index  and  FJ  and  F^  are  both 
evaluated  from  Equations  (4.1)  and  (4.2),  respectively. 

The  solution  is  obtained  by  iteration  since  the  values  of  Fxj 
depend  on  F'  T.  This  iteration  is  local  because  closure  is  obtained  at 
each  integration  point  before  proceeding  down  the  body  to  the  next 
integration  point.  Convergence  is  based  on  changes  of  less  than 
0.1  percent  in  both  the  heat  and  momentum  transfer  coefficients  between 
successive  Iterations.  If  the  iteration  fails  to  converge  In  30  tries,  a 
local  explicit  solution  is  obtained  by  setting  F*  T  «  F'  T  and 
subsequently  reevaluating  F'  .  based  on  the  resulting  value  of  Fv  T. 

A  |  1  X  y  l 

The  algorithm  then  proceeds  to  the  next  integration  point. 
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4.3  GEOMETRIC  ANALYSES 


It  i*  important  to  describe  the  geometry  with  high  resolution 
because  the  boundary  layer  integration  procedure  often  needs  small  step 
sizes  and  the  surface  metrics  and  normals  must  be  known  at  each 
integration  step.  This  surface  fitting  is  accomplished  in  3DMEIT  by  a 
parametric  interpolation  method  which  is  described  below  (Reference  26). 

The  vehicle  surface  is  separated  into  a  number  of  basic  elements, 
referred  to  as  patches.  Schematically,  patches  are  defined  as  having  four 
sides,  but  in  some  cases  one  of  the  sides  will  be  at  a  single  point.  The 
parametric  representation  of  a  patch  is  given  by: 

3  3 

r  =  [x(u,  w),  y(u,  w),  z(u,  w)]  -  ^  ]T  t..  u1'  w 1  (4.108) 

1-0  J-0 

where  0  <  u  <  1 

0  <  w  <  1 

Equation  (4.108)  transforms  the  surface  of  the  patch  into  a  unit  square  in 
(u  -  w)  space.  The  advantage  of  the  parametric  interpolation  scheme  is 
that,  once  identified,  the  patches  may  be  treated  independently,  and 
therefore,  only  a  single  patch  is  treated  in  detail  in  the  following 
discussion. 
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The  physical  identity  of  the  48  coefficients,  a^j,  which  must  be 
found  for  each  patch,  can  be  described  by  rewriting  Equation  (4.108)  in  the 
following  matrix  form: 

r  «  U  M  B  MT  WT  (4.109) 


where  the  superscript,  T,  refers  to  the  matrix  transpose,  and  where  U  and 
W  are  the  vectors  (u  ,  u  ,  u,  1)  and  (w  ,  w  ,  w,  1), 
respectively.  In  this  case,  the  matrix,  M  is  given  by: 

2  -2  1 

-3  3-2 

M  «s 

0  0  1 

1  0  0 

The  three  matrices,  B,  describe  the  input  quantities,  which  are  specified 
at  the  corner  points  of  the  patch. 


1 

-1 

0 

0 


(4.110) 


r{0, 

0) 

r(0,  1) 

I;'0-  » 

r(l, 

0) 

r(l,  1) 

I'1-  » 

2^ 

2- 

£(0 
3IJ ' 

,  0) 

(0,  0) 

’uSW(0-  » 

-fc. 

2- 

2- 

Kd 

3U  ‘ 

,  0) 

°> 

8  r  (1  0) 
3U3W'1’  ' 

UL(j  i) 

311 3W'  ’  ' 

(4.111) 


Calculating  the  elements  of  B  for  input  can  be  extremely  difficult.  The 
first  derivatives  would  require  that  the  parametric  form  of  the  curve  is 
explicitly  known  and  cross  derivatives  would  be  even  more  difficult  to 
calculate.  Therefore,  the  3DMEIT  code  uses  another  method  which  avoids 
inputting  derivatives. 
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The  48  coefficients  a^j  in  Equation  (4,108) are  computed  in  the 
code  using  the  input  locations  of  16  points  on  each  patch,  which  are 
usually  much  easier  to  generate  than  derivatives.  The  points  are  to  be 
uniformly  distributed  over  the  patch  surface  as  shown  in  Figure  4.3,  so 
the  input  points  will  map  into  a  uniform  grid  in  the  unit  square  in 
(u  -  w)  space.  Equation (4.109)  can  be  rewritten  in  the  following  manner: 


r  *  U  N  P  NT  WT 


(4.112) 


wnere 


and 


r(0,  0) 

r(0. 

1/3) 

r(0,  2/3) 

r(0,  1) 

f(l/3,  0) 

r( 1/3,  1/3) 

r(l/3,  2/3) 

r( 1/3,  1) 

(4.113) 

r ( 2/ 3 ,  0) 

r  (2/3 

,  1/3) 

r(2/3,  2/3) 

r(2/3,  1) 

r(l,  0) 

r(l, 

1/3) 

r(l, 

2/3) 

?(1,  1) 

-9/2 

27/2 

-27/2 

9/2 

9 

-45/2 

18 

-9/2- 

(4.114) 

-11/2 

9 

-9/2 

1 

1 

0 

0 

0 

Once  the  network  of  patch  points  has  been  developed,  the  code  gives 
designation  to  the  curves  by  requiring  that  the  vector  product  ar/au  x  ar/aw 

i 

produces  an  outward  surface  normal. 

Where  the  B  matrix  input  scheme  guarantees  continuity  in  slope  from 
one  patch  to  a  connecting  patch,  the  P  matrix  cannot  guarantee  this.  If 
slope  continuity  becomes  noticeably  absent,  this  is  remedied  within  the 
30MEIT  code  by  matching  the  unit  vectors  while  retaining  the  original 
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r 

s  ■ 


I 


9 


Jh 

magnitudes.  The  B  matrix  for  the  patch  Is  obtained  from  the  P  matrix  by 
the  following  op  mi: 


M*"1  N  P  NT 


(M"1) 


(4.115) 


If  the  slopes  in  the  u-direction  at  the  adjacent  corners  of  two  patches  do 
not  match,  continuity  is  assured  by  calculating  the  unit  vector  t  of  the 
slope  at  the  corner  of  the  first  patch  and  recalculating  the  value  of  the 
new  ar/au  for  the  second  patch.  The  new  value  of  ar/au  inserted  in  the  B 
matrix  for  the  second  patch  is  given  by: 


(4.116) 


The  new  B  matrix  is  then  used  to  determine  the  coefficients  a^  for  the 
patch.  Thus,  a  procedure  for  generating  a  parametric  interpolation 
surface  from  16  uniformly  spaced  points  has  been  established. 

This  procedure  allows  slope  discontinuities  to  exist  in  the  input  data 
that  define  the  corners  of  the  patches.  However,  the  internal  calculations 
described  above  remove  these  slope  discontinuities  and  facilitate  reliable 
boundary  layer  calculations.  This  is  necessary  because  boundary  layer  theory 
is  not  generally  applicable  to  bodies  with  slope  discontinuities. 

The  3DME1T  computer  code  limits  the  number  of  patches  on  a  vehicle 
to  60.  A  vehicle  is  set  up  in  segments  along  the  w  direction.  The 
maximum  number  of  segments  is  15.  For  each  segment  a  maximum  of  10 
patches  can  be  put  around  the  body,  but  the  overall  limit  of  60  patches 
must  be  observed.  Therefore,  when  15  segments  have  been  used  to  describe 
the  body,  only  four  patches  can  be  used  around  the  body  for  each  segment. 


W 


w 


» 
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If  the  vehicle  is  a  body  of  revolution,  four  points  along  a  meridional 
curve  can  be  input  for  each  segment  and  they  will  be  revolved  about  the 
axis  to  generate  the  required  number  of  patches.  Further  automation  of 
patch  generation  is  available  for  spherical  sections  and  certain 
maneuvering  vehicle  shapes  (see  Section  4.4, Input  Table  3). 

4.4  USER  INSTRUCTIONS 


This  user  oriented  section  describes  3DMEIT  inputs  and  outputs.  A 
complete  set  of  user  Instructions  is  presented  in  Section  4.4.1.  The 
3DME IT  output  format,  including  variable  descriptions,  is  presented  in 

Section  4.4.2.  Section  4.4.3  gives  a  sample  problem  and  selected  listings  of 
Input  and  output. 

4.4.1  INPUT  INSTRUCTIONS 

The  details  of  the  input  are  described  below.  The  basic  input  for 
a  run  consists  of: 

a  One  title  card 
•  Eight  input  tables 

Not  all  eight  of  the  input  tables  are  required  for  every  run.  Each  table 
is  preceded  by  a  single  card  containing  the  identifying  table  number. 

The  following  sections  describe  the  title  card,  and  the  eight  input 

tables. 

Title  Information 

The  first  card  contains  title  information  in  columns  1  through  80. 

The  contents  of  the  card  are  printed  at  the  beginning  of  the  boundary 
layer  solution. 
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The  following  table  describes  the  format  for  the  first  card  of  the 


data  deck. 

Card 

No. 

Columns 

Format 

Description 

Units 

1 

1-80 

8A10 

BCD  — 

Title  card  information 

(-) 

INPUT  ' 

TABLE  1. 

GENERAL  PROGRAM  CONSTANTS 

This  table  suppl 

ies  the  code  with  program  flags  that  indicate 

which 

options  to  be  subsequently  read 

and  other  computation  information. 

Card 

No. 

Columns 

Format 

Description 

Units 

1 

1-2 

12 

Enter  01  {table  number) 

(-) 

2 

1-5 

15 

NF (12) 
0 

1 

Plot  output  flag 

No  plots 

Read  Table  8  for  plot  instructions 

(-> 

6-10 

15 

NF (13) 
0 

1 

2 

Output  print  flag 

Summary  table  only 

Boundary  layer  output  for 
streamline  points 

Boundary  layer  output  for 
Integration  points 

(-) 

11-15 

15 

IRG 

0 

1 

Property  flag 

Ideal  gas  properties 

Real  gas  properties 

(-) 

16-20 

15 

NF  ( 5 ) 

1 

f  2 

Transition  criteria  flag 

Laminar  flow  (Table  04  not 
required) 

Critical  momentum  thickness 

(-) 

Table  04 
required 


^  4 

r  5 


Table  not 
required 


Reynolds  number  versus  edge  Mach 
number 

Critical  stream  length  Reynolds 
versus  edge  Mach  number 
Axial  distance  along  each  streamline 
Anderson  nose  criterion 


Re. 


&  I) 


0.7 


255  for  onset 
215  for  location 


*  "  VTe 

LORN  cone  criterion 
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INPUT  TABLE  1.  GENERAL  PROGRAM  CONSTANTS  (Continued) 


Card 

No.  Columns  Format  Description  Units 

6  Anderson  nose  criterion 

/ ^  0 . 7  255  for  onset 

^ee  \e  iji/  “  215  for  location 

=  0.  IB '  +  (0.9  +  0.1ICARB) 

(1  +  0.25B')(pe/pw) 

LORN  cone  criterion 

7  Fully  turbulent 

Positive  NREYCR  —  transitional 
heating 

Negative  NREYCR  —  no  transitional 
heating  (abrupt  transition) 

21-25  15  NF(6)  Carbon  flag  in  PANT  transition  (-) 

criterion  (NF(5)  -  6) 

0  Uses  0.9  factor  in  ^  calculation 
(recommended  for  carbon  materials) 

1  Uses  1.0  factor  in  ^  calculation 

26-30  15  NF (35 )  u-integratl on  mesh  density  (-) 

for  aerodynamic  coefficient 
calculations 

0  No  aerodynamic  coefficient 
calculation 

1  4  point  mesh  density 

2  8  point  mesh  density 

3  16  point  mesh  density 

31-35  15  NF ( 36 )  w-i ntegration  mesh  density  (-) 

for  aerodynamic  coefficient 
calculations 

0  No  aerodynamic  coefficient  calculation 

1  4  point  mesh  density 

2  8  point  mesh  density 

3  16  point  mesh  density 


57. 


INPUT  TABLE  I.  GENERAL  PROGRAM  CONSTANTS  (Continued) 


Card 


No. 

Columns 

Format 

Description 

Units 

2 

(cont.) 

36-40 

15 

IENTRP 

1 

2 

3 

4 

Entropy  swallowing  flag 

Normal  shock  entropy 

Inviscid  wall  entropy 

Streamtube  entrainment  option 

Local  mass  balance  option 

(-) 

41-45 

15 

I TWFLG 
-1 

0 

1 

Wall  temperature  flag 

Wall  temperature  calculated 

Constant  wall  temperature 
input  on  next  card 

Wall  temperature  distribution 
input  in  Table  6 

(-) 

46-50 

15 

IBPFLG 

-1 

0 

1 

Blowing  flag 

B'  calculated 

B'  -  0 

B*  distribution  input  in 

Table  6 

(-) 

51-55 

15 

JBLM 

0 

1 

Heat  transfer  blowing  correction 
flag 

Use  GE  option 

Use  Aerotherm  option 

(-) 

3 

1-10 

F10. 3 

SCALF 

1.0 

0.08333 

Scale  factor  for  input  configuration 

Input  in  feet 

Input  in  Inches 

11-20 

F10. 3 

TW 

Wall  temperature 

(°R) 

4 

This  card  is  necessary  only  if  NF ( 35 )  ^  0 
and  NF ( 36)  +  0 

1-10 

F10. 4 

DATAR(l)  Reference  area  for  aerodynamic 
coefficient  calculations 

(input2 

units) 

11-20 

F10.4 

DATAR{2)  Reference  length  for  aerodynamic 
coefficient  calculations 

(Input 

units) 

21-30 

F10. 4 

QATAR { 3 ) 

i  x-coordinate  of  the  center  of 
gravity 

( input 
units) 
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INPUT  TABLE  1.  GENERAL  PROGRAM  CONSTANTS  (Concluded) 


Columns  Format 


Description 


31-40  F10. 4  OATAR ( 4 )  y-coordlnate  of  the  center  of 

gravity 

41-50  FlO. 4  DATAR ( 5 )  z-coordinate  of  the  center  of 

gravity 


Units 

( i nput 
units) 

(input 

units) 
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INPUT  TABLE  2.  ENVIRONMENT  TABLE 


This  table  inputs  environment  conditions  and  the  streamlines  to  be 


used  for  the  boundary  layer  solution. 
Card 


No. 

Columns 

Format 

Data 

Units 

1 

1-2 

12 

Enter  02  (table  number) 

(-) 

2 

1-10 

F10. 2 

T1 

Freestream  temperature 

(°R) 

11-20 

F10. 2 

PI 

Freest.ream  pressure 

(atm) 

21-30 

F10.2 

VI 

Velocity 

(ft/s) 

31-40 

F10.2 

ALPHAB 

Angle  of  attack 

(deg) 

41-50 

no 

NSTREM 

Number  of  streamlines  (<  50) 

(-) 

3 

1-80 

8F10.2 

SANGLE 

(N),  N  -  1,  NSTREM,  Circum- 

(deg) 

ferenial  angles  of  the  stream¬ 
lines  at  the  sphere-cone 
tangency  point  (windward 
streamline  ■  180°) 
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INPUT 

TABLE  3. 

INITIAL  CONFIGURATION 

Card 

No. 

Columns 

Format 

Description 

Units 

1 

1-2 

12 

Enter  03  (table  number) 

(-) 

2 

1-10 

F10. 3 

RN 

Initial  nose  radius 

(in.) 

or 

(ft) 

11-20 

F10. 3 

THETC 

Initial  cone  half  angle 

(deg) 

21-30 

F10. 3 

BOA 

Height/span  ratio 

1.0  for  circular  cross  section 
t  1,0  for  elliptical  cross  section 

31-40 

F10. 3 

ZN 

Distance  from  apex  of  afterbody 
to  nosetip 

(in.) 

or 

(ft) 

41-45 

15 

NUMNOS 

Number  of  segments  on  the 
spherical  section  of  the  nosetip 

(-) 

46-50 

15 

NOSFLG 

Nosetip  generation  flag 

0  Each  segment  on  nosetip  must 
be  individually  input 

1  Code  automatically  sets  up 
NUMNOS  spherical  segments  on 
the  nose 

(-) 

3 

1-32 

1612 

NC  array 

(-) 

NC(1) 

Total  number  of  segments  (<15) 

NC(I), 

I  »  2,  NC < 1 )  input  mode  for  each 
segment  1  <,  NC  <  5 

If  NOSFLG  b  1,  values  must  be  input  for 
each  segment  but  the  code  will  ignore 
those  for  segments  2  through  NUMNOS 

4 

1-30 

1512 

NPATCH(I)  1=1,  15  number  uf  patches 

(-) 

around  the  body  on  segment  I 
(user  should  input  a  constant 
number  at  this  stage  in  code 
development) 


61. 


INPUT  TABLE  3.  INITIAL  CONFIGURATION  (Continued) 


The  optional  input  modes  for  the  segments  are  described  below.  The 
data  for  each  segment  must  be  ordered  to  match  the  mode  selection 
arrangement  specified  in  the  NC  array.  The  modes  can  be  separated  into 
three  options.  For  Option  1  (NC  <_  3),  a  single  meridian  curve  r(z)  is 
described  and  the  code  automatically  generates  NPATCH  patches  to  cover 
180°  circumferential ly  as  shown  in  Figure  4.4.  Option  2  is  used  to 
Input  completely  arbitrary  patches  as  also  shown  in  Figure  4.4,  Option  3 
is  used  to  describe  the  afterbody  of  a  vehicle  with  flat  sections  and  cut 
sections  as  illustrated  In  Figure  4.5.  The  input  formats  for  each  of 
these  options  are  as  follows: 


Input  Table  3(a).  Option  1 


NC ( I )  Format 


Description 


1  4F7.3  Conical  frustum,  input  zj,  rj, 

22*  r2 


r 


2 


3 


SF  7.3 


8F7.3 


Circular  arc  input  zj,  rj,  z2, 
r2«  02 


If  NOSFLG  =*  0,  each  segment 
cannot  turn  more  than  45° 


Arbitrary  curve  input 

21.  rl.  22,  r2,  z3,  r3,  Z4,  r4 

Note:  an  attempt  should  be 
made  to  distribute  evenly 
the  four  points  along  the 
curve  segment 
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INPUT  TABLE  3.  INITIAL  CONFIGURATION  (Continued) 
Table  3(b).  Option  2:  NC(I)  =  4 
Each  patch  requires  six  data  cards  which  are  ordered  as  follows: 
Card 


No._ 

Format 

Description 

1 

8F9.5 

x(0,  0),  x(0,  1/3),  x(U,  2/3),  x(0,  1),  x(l/3,  0), 
x(l/3,  1/3),  x( 1/3,  2/3),  x( 1/3,  1) 

2 

8F9.5 

x (2/3,  0),  x (2/ 3,  1/3),  x ( 2/ 3 ,  2/3),  x(2/3,  1), 
x(l,  0),  x(l,  1/3),  x(l,  2/3),  x(l,  1) 

3 

8F9.5 

y(0,0) ,  y(0,  1/3),  y(0,  2/3),  y(0,  1),  y(l/3,  0), 
y ( 1 / 3,  1/3),  y ( 1/ 3,  2/3),  y( 1/3,  1) 

4 

8F9.5 

y (2/3,  0),  y (2/ 3,  1/3),  y(2/3,  2/3),  y(2/3,  1), 
y(l,  0),  y(l,  1/3),  y ( 1 ,  2/3),  y(l,  1) 

5 

8F9.5 

z(0,0),  z(0,  1/3),  z(0,  2/3),  z(0,  1),  2(1/3,  0), 
z( 1/3,  1/3),  2(1/3,  2/3),  z(l/3,  1) 

6 

8F9.5 

2(2/3,  0),  z(2/3f  1/3),  2(2/3,  2/3),  z(2/3,  1), 
2(1,  0),  2(1,  1/3),  2(1,  2/3),  2(1,  1) 

Input  Table  3(a).  0^: :on  3:  NC(I)  =  5 


Columns 

Format 

Description 

Units 

1-7 

E7.3 

THETAC 

Cone  half  angle 

(deg) 

8-14 

E7.3 

VL 

Virtual  length  of  the 
measured  from  the  apex 
afterbody  cone 

vehicle 
of  the 

(in.) 

or 

(ft) 

15-21 

E7.3 

WCL 

Distance  from  the  end 
to  the  wi ndwr; d  cut 

of  the  body 

(in.) 

or 

(ft) 

22-28 

E7.3 

WF 

Width  of  the  flat 

(in.) 

or 

(ft) 
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INPUT  TABLE  3.  INITIAL  CONFIGURATION  (Continued) 


Columns 

Format 

Description 

Units 

29-35 

E7.3 

THETWS 

Angle  of  the  windward  cut 

(deg) 

36-42 

E7.3 

CLL 

Distance  from  the  end  of  the  body 
to  the  leeward  cut 

(in.) 

or 

(ft) 

43-49 

E7.3 

TilETSL 

Angle  of  the  leeward  cut 

(deg) 

50-56 

E7.3 

RS 

Radius  of  curvature  for  the  rounded 
corner  of  the  base 

(in.) 

or 

(ft) 

57-63 

E7.3 

BASXP 

Exponent  describing  the  cross 
sectional  shape  of  the  base 
=»  2.0  for  circular 
^  2.0  for  "super  circle" 

(-) 

64-65 

— 

Not  used 

66-70 

15 

K START 

Patch  number  where  the  interpolation 
between  an  exponent  of  2  and  BASXP 
begins 

(-) 

71-75 

15 

IFLAT 

Flat  flag 

0  No  flat  surface 

1  Flat  surface  on  windward  side 
determined  by  WF 

(-) 

The 

next  group  of  cards 

control  the  patch  setup  to  be  used  to  print 

the  boundary  layer  and 

surface 

energy  balance  results.  This  option  is 

for 

future  use 

of  3DMEIT  as  part  of 

a  procedure  for  predicting  shape  change.  It 

spline  fits 

a  group  of  patches 

in  a  least  squares  sense  and  then  redefines  the 

coordina tes 

of  the  body  segment  patches.  It  is  not  operative  in  the  present 

version  of 

3DMEIT,  but 

users  are  advised  to  enter  l's  in  all  fields  to 

prevent 

execution  errors  in  the  input  routines. 

Card 

No. 

Columns 

Format 

Description 

Units 

4  +  NC( 1) 

+  1 

1-30 

1512 

NDISC(I)  I  =  1,  15  The  number  of  w 
stations  per  spline  group 
to  be  used  in  the  recession 

(-) 

analysis 


INPUT  TABLE  3.  INITIAL  CONFIGURATION  (Concluded) 


No. 

Columns 

Format 

31-60 

1512 

61-62 

12 

4  +  NC ( 1 ) 

1-3 

13 

+  2 

4-6 

13 

4  +  NC ( 1 ) 

1-80 

8F10.4 

+  3 

4  +  NC( 1 ) 

1-80 

8F10.4 

+  4  »  n 

n  +  1 

1-80 

8F10.4 

n  +  2  »  m 

1-80 

8F10.4 

m.  +  1 

1-30 

3011 

m  +  2  »  i 

Description 

NDISC(I)  I  =  16,  30  The  number  of 
body  segments  per  spline 
group.  Up  to  15. 

Note:  it  is  suggested  that  the 

entire  nosetip  (spherical  section) 

be  included  in  one  group  and  the 

afterbody  in  another 

ND1SC(31)  Number  of  groups 

describing  the  body 

NTBZ  Number  of  w  locations 

where  boundary  layer  and 
surface  energy  balance 
results  are  saved  (30 
maximum) 

NTBU  Number  of  u  locations 

where  boundary  layer  and 
surface  energy  balance 
results  are  saved  (10 
maximum) 

TBZ(I)  I  *  1,  NTBZ  w  locations 

where  boundary  layer  and 
surface  energy  balance 
results  are  saved 

Continue  TBZ  array 

TBU(I)  I  =  1,  NTBU  u  locations 

where  boundary  layer  and 
surface  energy  balance 
results  are  saved 

Continue  TBU  array  input 

I MAT ( I ,  J)  1=1,  NTBZ  Surface 

material  flag  for  shape 
change  points 

Continue  for  J  =  2,  NTBU 


Units 

(-) 
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rl 


INPUT  TABLE  4.  TRANSITION  TABLE 


:rt 

Enter  this 

table  only  when  NF { 5 )  -  2,  3  or  4.  Maximum  of  50  entries 

are 

jr 

allowed  in 

this  table 

in 

1. 

Input  for  Re.versus  M„,  NF(5)  =.  2 

V  c 

i‘. 

Card 

§ 

No. 

Columns 

Format 

Description 

Units 

1 

1-  2 

12 

Enter  04  (table  number) 

(-) 

Nfc 

i  ■■■ 

1 

2 

1-  2 

12 

Nonzero  entry  for  the  last  card 
in  the  table 

(-) 

i 

lk‘  •' 

s  , 

1 

3-14 

E12.5 

AM  Local  edge  Mach  No.  (Me) 

(-) 

f,'; 

[.' 

15-26 

El. 5 

REM  Critical  momentum  thickness 

(-) 

.* 

ir 

Reynolds  number  (Re0) 

k 

3 

Same  as  Card  2  for 

increasing  Me 

(etc.) 

V: 

2. 

Input  for  Reg  versus  M0,  NF ( 5 )  *  3 

*  *.  1 

Card 

■ 

a 

No. 

Columns 

Format 

Description 

Units 

»i 

1 

1-  2 

12 

Enter  04  (table  number) 

(-> 

4 

2 

1-  2 

12 

Nonzero  entry  for  the  last  card 
in  the  table 

(-> 

\ 

\ 

.  *i 

’ » 

3-14 

E12.5 

AM  Local  edge  Mach  No.  (Me) 

(-) 

•l 

**  i 

;■  ■■ 

15-26 

El. 5 

REM  Critical  stream-length 

(-) 

■  1 
•  *< 

!-;■ 

Reynolds  number  (Res) 

■  i 

:  «' 

r- 

3 

Same  as  Card  2  for 

increasing  Me 

*  ’ 

(etc.) 
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3.  Input  for  Transition  Location,  NF ( 5 )  *  4 


Columns 

Format 

Description 

1-  2 

12 

Enter  04  (table  number) 

1-  2 

12 

Nonzero  entry  for  the  last  card 
in  the  table 

3-14 

£12.5 

AM  Axial  location  of  transition 
along  each  streamline  specified 
in  Table  2 

INPUT  TABLE  5.  MATERIAL  PROPERTIES 


Maximum  of  three  material  inputs  are  allowed  in  this  table.  Each 
material  input  is  treated  according  to  the  following  format. 


Card 

No. 

Columns 

Format 

Description 

Units 

1 

1-  2 

12 

Enter  05  (table  number) 

(-) 

2 

1-  5 

15 

MAT 

Material  index  for  following 
table 

(-) 

6-10 

15 

UROUGH 

Laminar  heating  augmentation 
flag 

0  No  augmentation 

1  Transition  proximity 
augmentation 

(-) 

3 

1-  2 

Blank 

3-14 

E12.5 

RUFL 

Intrinsic  roughness  height; 
nonuniform  roughness  may  be 
input  in  Table  6 

(mil) 

15-26 

E12.5 

RUFMAX 

Maximum  turbulent 
roughness  height  (k.  ) 

lmax 

(mil) 

27-38 

E12.5 

RUF1 

Constant  ki  (mi l-psi^*77) 

>0  White-Grabow  scallop 
law  is  used: 

k*  »  k,  p  ~0,77  with  k.  -  RUFMAX 

t  1  e  t 

=0  Constant  turbulent  roughness  of 

RUFMAX  is  used 

39-50 

Ei2.5 

EMISS 

Emissivity 

(-) 

51-62 

E12.5 

CMH 

Ratio  of  mass  to  heat 
transfer  coefficients 

(-) 
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INPUT  TABLE  5.  MATERIAL  PROPERTIES  (Concluded) 


Parameters 

Columns 
1-  2 
3-14 
15-26 

27-38 

39-50 


n  Blowing  Correction  to  Transfer  Coefficients 

Def aul t 

Format  Description  Values 


12 

NC 

Flag,  zero 

E12.5 

BLS 

Laminar  shear  parameter 

0.5 

E12.5 

BLH 

Laminar  heating 
parameter 

0.5 

E12.5 

BTS 

Turbulent  shear 
parameter 

0.35 

E12.5 

BTH 

Turbulent  heating 
parameter 

0.35 

If  zero  is  entered,  default  values 
are  used. 


INPUT  TABLE  6.  SURFACE  DATA 


The  3DMEIT  code  Includes  approximate  techniques  to  predict  surface 
temperature  and  blowing  rate  distributions  of  ablating  surfaces.  The 
user,  however,  may  bypass  these  computations  and  enter  the  desired  surface 
temperature  or  blowing  rate  distributions  via  this  table.  A  nonuniform 
surface  roughness  distribution  may  also  be  input  in  this  table.  Entries 
for  each  variable  are  described  in  the  following  subsections.  Any 
combination  of  one  or  three  subsections  can  be  input. 

Card 


No. 

Columns 

Format 

Description 

Units 

1 

1-2 

12 

Enter  06  (table  number) 

(-) 

2 

1-5 

15 

ITABL  Enter  l  (temperature  table) 

(-) 

3 

1-80 

9F10.5 

Surface  temperature  for  NTBZ  w  points 

(°R) 

4 

1-80 

8F10.5 

Same  as  above  for  each  of  NTBU  u  values 

(°R) 

m 

1-5 

15 

ITABL  Enter  2  (blowing  rate  table) 

m+1 

1-80 

8F10.5 

Blowing  rate  (B1)  »  Pwvw/Peuecm 
for  each  NTBZ  w  values 

(-) 

m+2 

1-80 

8F10.5 

Same  as  above  for  each  NTBU  u  value 

(-) 

n 

1-5 

15 

ITABL  Enter  3  (surface  roughness) 

n+1 

1-80 

8F10.5 

Surface  roughness  for  NTBZ  w  points 

(mil) 

n+2 

1-80 

8F10.5 

Same  as  above  for  each  of  NTBW  w  values 

(mil) 

Last 

card 

1-5 

15 

ITABL  <  0  signifies  end  of  Table  6 

(-) 

INPUT  TABLE  7.  INVISCID  DATA 


The  3DMEIT  code  reads  inviscid  data  from  tape  23  and  stores  the 
needed  edge  conditions  and  integrated  mass  fluxes  on  tape  22.  If  a 
calculation  is  to  be  done  using  the  same  inviscid  data  as  a  previous 
calculation,  one  may  use  the  stored  values  on  tape  22  and  avoid  the  cost 


of  integrating 

the  velocity  fields 

Card 

No. 

Columns 

Format 

1 

1-  5 

15  IRSTRT 

For 

IRSTRT  *  0, 

the  format  on  tape 

Card 

No. 

Columns 

Format 

1 

1-13 

E13.6 

AMINF 

14-26 

E 13. 6 

ALPHA 

27-29 

13 

MDAT 

30-32 

13 

NDAT 

2 

1-13 

E13.6 

7.  DAT 

3  > 

3  +•  MDAT 

1-13 

E13.6 

PHDAT 

14-26 

E13.6 

RB 

and  calculating  the  metrics. 

Description  Units 

Inviscid  data  input  flag  (-) 

0  Read  inviscid  data  from 
tape  23 

1  Read  previously  defined  mass 
fluxes  and  metrics  from 
tape  22 

23  is  given  in  the  following  table: 


Description 

Units 

Freestream  Mach  number 

(-) 

Angle  of  attack 

(deg) 

Number  of  steamlines 

(-) 

Number  of  points  through 
the  shock  layer 

(-) 

Axial  location 

(in.) 

Circumferential  angle 

(rad) 

Body  radius 

(in.) 

arg/aZ 

(-) 

'  -• 


«r 


V 


27-39 


E13.6  RBZ 
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INPUT  TABLE  7.  INVISCID  DATA  (Concluded) 


Columns 

Format 

Description 

Units 

♦  MDAT  (cont.) 

40-52 

E13.6 

RBTH 

3rg/BQ 

(in.) 

53-65 

E13.6 

OMEGA 

tan“l(u^/uz)  on  the 
surface 

(rad) 

66-78 

E13.6 

SF 

Stream  function 

(-) 

+  (NDAT  * 

MDAT) 

1-13 

E13.6 

PF 

Flowfield  pressure 
p/p^  IRG  *  0 
p  IRG-1 

(-) 

(lbf/ftZ) 

14-26 

E13.6 

T 

Temperature 

T/T  IRG  -  0 

T  "  IRG  a  1 

(-) 

(°R) 

27-39 

E13.6 

VR 

Radial  velocity  component 

u„  ,/p*  /p  IRG  -  0 
r  v  oo  oq 

ur  IRG  «  1 

(-) 

(ft/s) 

40-52 

E13.6 

VP 

Circumferential  velocity 
component 

/ p  /p  IRG  -  0 

uy  00  *  IRG-1 

(-) 

(ft/s) 

53-65 

E13.6 

VZ 

Axial  velocity  component 

uzJp  /pm  IRG  -  0 
uz  *  00  IRG-1 

(-) 

(ft/s) 

66-78 

E13.6 

YR 

Distance  from  body 

(in.) 

Card  sets  2,  3,  4  are  repeated  for  each  axial  station. 
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INPUT  TABLE  8. 


PLOT  DATA 


Card 

No.  Format  Description 


Units 


1 


12  Enter  08  (table  number) 


(-) 


F7.3 

YL 

Minimum  value  of 

F7.3 

YR 

Maximum  value  of 

F7.3 

XB 

Minimum  value  of 

F7.3 

XT 

Maximum  value  of 

F7.3 

ZN 

Minimum  value  of 

F7.3 

ZB 

Maximum  value  of 

y 

on 

the 

vehicle 

[ft 

or 

in. ) 

y 

on 

the 

vehicle 

tft 

or 

in. ) 

X 

on 

the 

vehicle 

[ft 

or 

in. ) 

x 

on 

the 

vehicle 

[ft 

or 

in. ) 

z 

on 

the 

vehicle 

[ft 

or 

in. ) 

z 

on 

the 

vehicle 

!ft 

or 

in. ) 

Note:  these  six  dimensions  should  form  a  box 
completely  enclosing  the  vehicle  to  aid  the 
system  In  scaling  the  plot 


12 

NU 

Number  of  u-increments  to  be  taken  < 

10 

(-) 

12 

NW 

Number  of  w-increments  to  be  taken  < 

10 

(-) 

12 

NPLOT 

Number  of  views  at  a  given  time 

<-) 

l' 0  draws  complete  grid 

12 

NI  -  l 

1  draws  only  patch  boundaries 

(-) 

^  2  draws  profi 1e 

F7.3 

XBETA 

View  orientation  angle  defined  In 
Figure  4.6 

(deg) 

F7.3 

XTHETA 

View  orientation  angle  defined  in 

(deg) 

Figure  4.6 


4  >  n 


Same  as  card  number  3  for  NPLOT  views 


End  of  Input 

The  end  of  the  input  deck  is  signaled  by  a  single  card  with  a  -1 
punched  in  columns  1  and  2. 
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4.4.2  OUTPUT  DESCRIPTION 

This  section  describes  the  3DMEIT  code  output.  There  are  two 
classes  of  output  from  the  3DMEIT  code: 

1.  Printout  of  the  input  data 

2.  Calculation  (results)  output 
These  are  described  in  the  following  sections. 

4.4.2. 1  Printout  of  Input  Data 

Program  output  begins  with  a  complete  printout  of  the  contents  of 
input  Tables  1  through  7.  Most  of  this  output  is  fully  labeled  and  is 
printed  exactly  as  input  by  the  user. 

4 . 4 . 2 . 2  Calculation  Output 

The  results  of  the  calculation  are  output  In  one  of  three  formats. 

All  three  output  formats  Include  up  to  five  separate  groups  of  Information, 

These  groups  are  (1)  environment  information,  (2)  fluid  properties, 

(3)  boundary  layer  results,  (4)  body  geometry  and  surface  quantities 
results,  and  (5)  aerodynamic  coefficients.  The  output  formats  differ  in  the 
degree  of  detail  of  the  boundary  layer  results.  Also,  the  aerodynamic 
coefficients  are  only  printed  if  NF (35)  and  NF ( 36)  are  nonzero.  The 
following  sections  describe  the  variables  printed  in  each  of  the  output 
groups. 

ENVIRONMENT  INFORMATION 

•  ENVIRONMENT  NUMBER 

•  DESCRIPTIVE  HEADING:  input  on  title  card 
■  •  MACH  NO:  current  freestream  Mach  number 

•  VELOCITY:  current  freestream  velocity  (ft/s) 

•  ANGLE  OF  ATTACK:  current  angle  of  attack  (deg) 

t  UNIT  RE  NO:  current  value  of  freestream  unit  Reynolds  number  (per  ft) 


•  'LUID  PROPERTIES 

The  following  results  are  printed  for  the  freestream  conditions, 
the  conditions  behind  the  normal  shock,  and  the  stagnation  point 
conditions: 

•  PRESSURE  (atm) 

•  TEMPERATURE  (°R) 

•  DENSITY  (lbin/ft3) 

•  ENTHALPY  (Btu/lbm) 

•  ENTROPY  (Btu/lbm  °R) 

•  ISENTROPIC  EXPONENT 
BOUNDARY  LAYER  RESULTS 

The  summary  output  ( NF ( 13 )  ■  0)  includes  only  the  first  of  the 
following  sots  of  information.  The  rest  of  the  tables  a. e  printed  for 
each  streamline  at  either  the  streamline  points  ( NF ( 13)  *  1)  or,  in  more 
detail,  at  the  boundary  layer  integration  points  (NF(13)  a  2). 

Boundary  Layer  Summary 

•  STREAMLINE  NO:  streamline  number 

•  BETA,  deg:  initial  direction  of  the  streamline  from  the 
stagnation  point 

i e  SONIC  POINT  LOCATION,  in.:  stream! ength  from  the  stagnation 

point  to  the  sonic  point,  and  the  Cartesian  coordinates 

•  TRANSITION  POINT  LOCATION,  in.:  streamlength  from  the 
stagnation  point  to  the  transition  ooint,  and  the  Cartesian 
coordi nates 

•  POINT  NO:  streamline  point  number  or  boundary  layer 
integration  point  number 
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•  STREAMLINE  LOCATION,  in.  and  deg:  z-coordinate  (measured  from 
initial  z-location  of  the  nosetip)  and  PHI  —  the  angular 


location  on  the  body  (measured  from  the  positive  x-axis  in  the 
cross  sectional  plane) 

•  METRIC  COEFFICIENT,  in.:  metric  coefficient  for  the  streamline 
coordinate  system 

•  BODY  ANGLE,  deg:  angle  that  a  line  tangent  to  the  surface  and 
perpendicular  to  the  normal  makes  with  the  wind  vector 

•  EDGE  PRESSURE,  atm:  pressure  at  the  edge  of  the  boundary  layer 

•  WALL  TEMPERATURE,  °R:  temperature  at  the  wall 

i  B-PRIME  THERMOCHEM:  dimensionless  blowing  parameter, 

B‘  =  (pv)  /M  due  to  ablation  only,  where  M  is  the  mass 

W 

transfer  coefficient 

•  RECOVERY  ENTHALPY,  Btu/lbm:  hr  £  hg  +  r  (ht  -  hg) 
where  ht  is  the  total  enthalpy,  hg  is  the  edge  enthalpy, 
and  r  is  the  recovery  factor 

2 

•  HEAT  TRANS  COEFFICIENT,  lbm/ft  -s:  H  *  PeugCh  - 

r  w 

where  C.  is  the  Stanton  number  and  q, ,  is  the  surface  heat  flux 
n  ^w 

Viscous  Flow  —  Wall  and  Boundary  Layer  Recovery  Properties 


• 

BODY  PT  NO  (J):  index  of  the  body  points 

* 

• 

INTEG  PT  NO  (I):  index  of  the  integration  point  for 

computed  parameters  are  printed 

which  the 

« 

t 

WALL  TEMPERATURE,  °R  (TW):  temperature  of  the  wall 

• 

■ 

WALL  ENTHALPY,  Btu/lbm  (HW):  enthalpy  of  the  flow  at 

temperature 

the  wall 

1  ■ 

• 

WALL  DENSITY,  lbm/ft3  (ROW):  density  of  the  flow  at 

temperature 
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the  wall 

'  V 

:  *  . 

►  - 

m 

i 

•1 


•  WALL  VISCOSITY,  lbm/ft-s  (VISW):  viscosity  at  the  wall 

•  RECOVERY  ENTHALPY,  Btu/lbm  (HR):  recovery  enthalpy, 

hr  “  he  *  r(ht  '  he} 

•  RECOVERY  FACTOR  (HECOV):  recovery  factor,  r 

•  SENSBL  CONV  HEAT  FLUX,  Btu/ft2-s 

•  CF/2:  skin  friction  coefficient  divided  by  two, 

2 

Cf /2  ■  t^/ PgUg 

•  B-PRIME  THERMOCHEM:  dimensionless  blowing  parameter, 

B'  =  (pv)  /M  due  to  ablation  only,  where  M  is  the  mass 

W 

transfer  coefficient 

•  RECOVERY  ENTHALPY,  Btu/lbm:  hr  =  hfi  +  r  (ht  -  hfi)  where  ht  is 
the  total  enthalpy,  hg  is  the  edqe  enthalpy,  and  r  is  the 
recovery  factor 

2  qw 

•  HEAT  TRANS  COEFFICIENT,  lbm/ft-s:  H  -  peugCh  - 

where  is  the  Stanton  number  and  qw  is  the  surface  heat  flux 
Viscous  Flow  —  Edqe  Properties  Table 

•  BODY  PT  NO  (J):  index  of  the  body  parts 

•  INTE6  PT  NO  (I):  index  of  the  integration  point  for  which  the 
computed  parameters  are  printed 

•  STREAM  LENGTH,  in.  (S):  stream  length  along  the  nosetip 
surface  from  the  stagnation  point  to  the  integration  points 

•  VELOCITY,  ft/s  (UE):  velocity  at  the  edge  of  the  boundary  layer 

•  MACH  NO  (HCAM) :  Mach  number  at  the  edge  of  the  boundary  layer 

•  ENTHALPY,  Btu/lbm  (HE):  enthalpy  at  the  edge  of  the  boundary 
layer 
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•  TEMPERATURE,  °R  (TE):  temperature  at  the  edge  of  the 
boundary  layer 

•  DENSITY,  lbm/ft^  (ROE):  density  at  the  edge  of  the  boundary 
1  ayer 

•  VISCOSITY,  lbm/ft-s  (VISE):  viscosity  at  the  edge  of  the 
boundary  layer 

0  UNIT  RE  NO,  1/ft  (URE):  unit  Reynolds  number  at  the  edge  of 
the  boundary  layer 

Viscous  Flow  —  Boundary  Layer  Solution  Table 

0  BODY  PT  NO  (J):  index  of  the  body  points 
0  INTEG  PT  NO  (I):  index  of  the  integration  point  for  which  the 
computed  parameters  are  printed 
0  STREAM  LENGTH,  in.  (S):  stream  length  along  the  nosetip 
surface  from  the  stagnation  point  to  the  integration  points 
0  MOMENTUM  THICKNESS,  mil  (THE):  momentum  thickness  e  of  the 
boundary  layer 


where  6  is  the  momentum  boundary  layer  thickness  and  y  is  the 
radial  distance  from  the  centerline.  For  a  more  general 
definition  of  e  (also  i  and  6*)  see  Reference  1). 

0  ENERGY  THICKNESS,  mil  (PHI):  energy  thickness  t  of  the 
boundary  layer 
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where  ht  is  the  total  enthalpy  and  is  thickness  of  the 
energy  boundary  layer 

•  SHAPE  FACTOR  (HSF):  boundary  layer  shape  factor,  H  »  a*/e 

•  MOM  THICK  RE  NO  (RETH):  Reynolds  number  based  on  the  momentum 
thickness 


•  HEAT  TRANS  COEFFICIENT,  lbm/ft2s  (RUCH):  heat  transfer 
coefficient,  H  «=  peueCh  «  qw/(hr  -  hw) 

where  C^  is  the  Stanton  number  and  qw  is  the  wall  heat  flux 

•  REYNOLDS  ANAL  FAC  (RAF):  Reynolds  analogy  factor 

Ch/(l/2)Cf, 

o 

where  C../2  •-  tu/(p0u0  )  is  the  friction  factor  and 

T  W  C  c 

a-  is  the  wall  shear 

W 

•  INTERMITTENCY  (ADML):  boundary  layer  intermittency  factor  f, 
where  0  <  f  <  1  for  flow  ranging  from  fully  laminar  to  fully 
turbulent  regime 

•  TRANSITION  PARAMETER:  Anderson  nosetip  transition  parameter. 


where  ^  =  0.1  B'  +  0.9  (1  +  0.25  B')  p e/pw 


R2. 


Viscous  Flow  —  Curved  Shock  and  Roughn ess  Effects  Table 

•  BODY  PT  NO  (0):  Index  of  the  body  points 

•  INTEG  PT  NO  (L):  index  of  the  integration  point  for  which  the 
computer  parameters  are  printed 

•  STREAM  LENGTH,  in.  (S):  stream  length  along  the  nosetip 
surface  from  the  stagnation  point  to  the  integration  points 

•  EDGE  ENTROPY,  Btu/lbm  -  °R  (ENTR):  entropy  at  the  edge  of 
the  boundary  layer 

•  EDGE  STREAMLINE  LOCATION  AT  SHOCK,  in.  (YBAR) :  radial 
coordinate  of  the  point  of  intersection  of  the  boundary  layer 
edge  streamline  with  shock  wave 

•  EDGE  MASS  FLUX  AUGMENTATION  (ROUE):  ratio  of  the  boundary 
layer  edge  mass  flux,  p  u  to  that  for  a  normal  shock 

W  C 

•  ROUGHNESS,  mil  (RUF):  surface  effective  sand  grain  roughness 
height 

•  HEAT  TRANSFER  AUGMENTATION  (RUFSMT) :  augmentation  of  the  heat 
transfer  coefficient  due  to: 

—  transition  proximity  for  the  laminar  portion,  and 

—  surface  roughness  effects  for  the  turbulent  portion  of  the 
boundary  layer 

BODY  GEOMETRY  AND  BOUNDARY  LAYER  RESULTS 

This  table  includes  information  for  all  of  the  u  and  w  locations 
where  the  information  is  saved  (TBU  and  TBZ  arrays,  see  Table  3,  input). 

•  STAGNATION  POINT  LOCATION:  Cartesian  coordinates  and  the 
normalized  surface  coordinates  of  the  stagnation  point 


•  MINIMUM  Z-COORD  LOCATION:  Cartesian  coordinates  and  the 
normalized  surface  coordinates  of  the  point  with  the  minimum 
z-coordinate 

•  SURFACE  COORDINATES:  non-normal i zed  u  and  w  locations  of  the 
points 

•  CARTESIAN  COORDINATES:  Cartesian  coordinates  of  the  points 

•  RECOVERY  ENTHALPY  (same  as  boundary  layer  output) 

•  SURFACE  PRESSURE  (same  as  edge  pressure  in  boundary  layer) 

•  HEAT  TRANS  COEFFICIENT  (same  as  boundary  layer  output) 

•  INDUCF.D  PRESSURE:  pressure  calculated  using  the  boundary 
displacement  thickness  (atm) 

•  SKIN  FRICTION:  skin  friction  coefficient,  Cf/2 

•  B-PRIME  THERMOCHEM  (same  as  boundary  layer  output) 

•  SURFACE  TEMPERATURE  (same  as  boundary  layer  output) 

AERODYNAMIC  COEFFICIENTS 

This  table  is  printed  if  the  aerodynamic  coefficients  are 
calculated  (NF (35)  ^  0  and  NF ( 36 )  ^  0).  The  output  Includes  the 
lift-to-drag  ratio,  the  center  of  pressure,  and  the  x,  y,  and  z  components 
of:  (1)  drag,  (2)  lift,  (3)  total  force,  (4)  inviscid  force,  (5)  viscous 
forces  (induced  pressure  and  shear),  (6)  total  moment,  (7)  inviscid 
moment,  and  (8)  viscous  moments. 

4.4.3  SAMPLE  PROBLEM 

This  section  presents  input  listings  and  selected  output  page 
listings  for  the  flight  test  case.  The  complete  output  is  shown  for  the 
windward  and  leeward  streamlines,  while  only  the  first  output  table  is 
shown  for  the  remaining  streamlines. 
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SECTION  5 


RESULTS 


The  two-dimensional  version  of  ME  IT  has  boon  verified  by  many 
comparisons  with  data  and  other  calculations  (References  6  told). 
Therefore,  testing  of  3DMEIT  has  been  concentrated  on  data  from  hi  coni  c 
configurations  at  angle-of-attack.  Three  runs  from  References  27  and  28 
were  selected  as  test  cases  for  3DME IT.  The  conditions  of  these  three 
test  cases  are  summarized  in  Table  5-1.  A  demonstration  calculation  was 
also  done  at  flight  conditions  on  a  biconi;  configuration  with  yaw 
stabilizers.  This  calculation  exercised  the  ablation  model  to  predict 
wall  temperature  and  blowing  parameters. 

Figure  5-1  compares  an  early  zero  angle-of-attack  heat  transfer 
prediction  with  the  experimental  data.  In  general,  this  prediction  is 
higher  than  the  data,  which  is  believed  to  be  due  to  the  ideal  gas 
properties  used  in  this  calculation.  The  local  calculated  stagnation 
temperature  was  1,982°R,  which  is  higher  than  the  measured  stilling 

i 

chamber  temperature  of  1,894°R.  The  measured  local  freestroam 
stagnation  temperature  was  1 , 7 9 6 0 R  which  indicates  some  losses  from  the 
stilling  chamber  conditions.  Therefore,  the  freestream  conditions  for  the 
calculations  were  modified  as  shown  in  Table  5-1  to  more  closely  match  the 
measured  local  stagnation  temperature.  This  resulted  in  much  better 
agreement  with  the  experimental  data,  as  shown  in  subsequent  comparisons. 
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Figures  5-2  and  5-3  show  the  laminar  calculation  of  5°  angle-of-attack. 
Again,  a  stagnation  temperature  of  1,998°R  gives  predictions  greater  than 
the  data.  However,  with  the  modified  stagnation  conditions,  the  predic¬ 
tions  are  quite  good,  except  on  the  leeward  side.  Figure  5-3  shows  that 
this  discrepancy  occurs  where  viscous  interactions  would  be  expected  to 
affect  the  inviscid  pressure. 

The  solution  for  turbulent  flow  at  a  =  10°  is  shown  in  Figures  5-4 
through  5-7.  Figure  5-4  illustrates  the  effect  of  entropy  swallowing 
and  shows  a  significant  increase  in  heating  on  the  windward  side 
compared  to  isentropic  flow.  The  predictions  with  entropy  swallowing 
are  in  good  agreement  with  the  experiment. 

Like  the  laminar  predictions,  those  for  turbulent  flow  are  poorer 
on  the  leeward  side.  Figures  5-6  and  5-7  compare  the  inviscid  pressure 
and  the  experimental  data.  Again,  the  inviscid  procedure  underpredicts 
the  pressure  on  the  forecone  and  overpredicts  the  pressure  on  the  aft 
cone  because  of  viscous  interaction  effects  on  the  leeward  side. 

The  fourth  case  is  a  demonstration  of  the  flight  capabilities  of 
3DMEIT.  This  case  used  real  gas  properties,  ablation  modeling  and  a 
complex  configuration  with  yaw  stabilizers  on  the  aft-cone.  The  patch 
geometry  used  in  the  calculation  Is  shown  in  Figure  5-8.  Figures  5-9 
and  5-10  give  the  heat  transfer  on  this  configuration. 
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Figure  5-3.  Circumferential  Heat  Transfer  Distributions,  Groups  96 
97,  136  to  158  Laminar,  M  «  10,  o  =  5° 
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Figure  5-6.  Inviscid  Surface  Pressure  Distribution 
Compared  with  HYTAC  Data,  M  =  8,  a  =  10° 
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Figure  5-7.  Inviscid  Surface  Pressure  Distributions 
Compared  with  HYTAC  Data,  M  *  8,  a  =  10° 


IBtu 


V  =  15,609  ft/sec 
M  =  16.2 

CD 

p  ■  0.075  atm 

roo 

T  =  389. 87°R 


a  =  6 


Alt  -  6,000  ft 


Axial  distance  (in.) 


Figure  5-9.  Predicted  Axial  Heating  Distributions  at  Flight  Conditions, 
10. A/6.0  Biconic  with  Yaw  Stabilisers,  M  =  16.2,  a  =  5° 
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SECTION  6 


CONCLUSIONS 

This  report  documents  the  initial  efforts  In  the  development  of 
a  coupled  Inviscld-boundary  layer  flow  field  solution  procedure  speci¬ 
fically  designed  to  treat  maneuvering  reentry  vehicles.  The  resultant 
code,  which  consists  of  the  3IS  inviscid  flow  field  code,  a  buffer  code, 
and  the  3DMEIT  Integral  boundary  layer  code,  can  treat  equilibrium  air 
thermodynamics,  laminar,  transitional,  or  turbulent  flows,  and  models 
the  effects  of  surface  roughness  and  mass  addition  from  ablating  heat- 
shields  on  the  boundary  layer  flow.  Approximate  models  are  provided  to 
automatically  determine  the  wall  temperature  and  mass  addition  rates  over 
the  entire  vehicle,  for  a  variety  of  heatshleld  materials. 

This  coupled  inviscld-boundary  layer  flow  field  procedure 
represents  an  extension  to  the  capabilities  of  such  techniques,  obtained 
by  making  maximum  use  of  the  detailed  inviscid  flow  field  information 
available.  Specifically,  this  approach  has  allowed: 

•  Direct  calculation  of  the  inviscid  surface  streamlines 

in  the  finite-difference  inviscid  afterbody  code  (3IS), 
and 

•  Accurate  definition  of  the  boundary  layer  edge  properties 

through  the  use  of  the  detailed  inviscid  flow  field 
profiles  near  the  wall,  rather  than  using  a  stream  tube 
mass  balance  technique  (which  is  difficult  to  apply  to 
a  three-dimensional  flow). 


The  flow  field  technique  resulting  from  this  effort  provides  a 
reliable  method  for  engineering  predictions  of  aerodynamic  loads  and  heat¬ 
ing  on  maneuvering  reentry  vehicles,  while  being  efficient  enough  (through 
the  use  of  an  integral  boundary  layer  procedure)  to  allow  its  routine  use 
iri  both  vehicle  design  and  analysis. 

It  is  recommended  that  future  efforts  be  conducted  to  1)  extend 
the  3DME IT  method  for  body  geometry  description  in  order  to  enable  calcu¬ 
lations  to  be  performed  at  sideslip  and  2)  allow  for  an  interface  between 
the  buffer  code  and  the  output  of  a  CM3DT  solution  in  order  to  treat 
non-spherical ,  asymmetric  nosetip  shapes. 
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